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profile.  The  results  are  applicable  to  rough  surfaces  which  do  not  have  rms  slopes  greater  than  a 
few  degrees. 

The  solution  is  applied  to  laser  scattering  from  aluminum  and  titanium  airplane  surfaces  where  the 
medium  is  the  clear  atmosphere.  Proflcorder  tracing  of  sample  rough  surfaces  reveals  independent 
"roughness^and  "wavlness"  components  normally  distributed  with  Gaussian  (isotropic)  auto¬ 
correlation.  The  deterministic  surface  model  is  a  generalized  quadratic  approximation  to  the 
surface  mean.  The  laser  profile  is  the  lowest-order  circular  laser  mode.  The  atmospheric  OTF  Is 
modeled  as  Gaussian. 

For  typical  parameters,  the  result  decouples  Into  a  specular  term  (depending  on  turbulence,  wave 
profile,  and  the  known  surface)  and  a  diffuse  component  (depending  on  the  rough  and  known  surfaces). 
The  deterministic  surface  Is  deformed  spherically  to  account  for  an  Incident  wave  of  slight  spherical 
character.  (A  spherical  wave  affects  only  the  coherent  scatter.)  From  a  measurement  of  scatter  at 
normal  incidence  the  surface  power  spectral  density  and  the  atmospheric  refractive  index  structure 
constant  may  be  determined  provided  the  surface  is  isotropic,  homogeneous ,  and  normally  distributed. 

Rough  surface  and  profile  aspects  of  the  theory  are  verified  by  nearfleld  measurements  of  scatter 
from  titanium  and  aluminum  when  there  is  no  turbulence  present.  Surface  power  spectral  densities 
found  empirically  from  scatter  at  normal  Incidence  agree  with  the  spectra  found  from  proflcorder 
tracing.  i 
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ABSTRACT 


A  generalized  investigation  of  wave  scattering  from  rough  surfaces 
yields  an  expression  for  mean  scattered  power  flux  which  is  the  spatial 
Fourier  transform  of  the  product  of  three  transfer  functions:  (a)  the  optical 
transfer  function  (OTF)  of  the  propagating  medium,  (b)  the  joint  characteristic 
function  of  the  surface  roughness,  and  (c)  a  coherence  involving  the  mean 
surface  and  the  wave  profile.  The  results  are  applicable  to  rough  surfaces 
which  have  rms  slopes  less  than  a  few  degrees. 

The  solution  is  applied  to  laser  scattering  from  aluminum  and  titan¬ 
ium  airplane  surfaces  where  the  medium  is  the  clear  atmosphere.  Proficorder 
tracing  of  sample  rough  surfaces  reveals  independent  "roughness"  and  "wavi¬ 
ness"  components  normally  distributed  with  Gaussian  (isotropic)  autocorrela¬ 
tion.  The  deterministic  surface  model  is  a  generalized  quadratic  approximation 
to  the  surface  mean.  The  laser  profile  is  the  lowest-order  circular  laser  mode. 
The  atmospheric  OTF  is  modeled  as  Gaussian. 

For  typical  parameters,  the  result  decouples  into  a  specular  term 
(depending  on  turbulence,  wave  profile,  and  the  known  surface)  and  a  diffuse 
component  (depending  on  the  rough  and  known  surfaces) .  The  deterministic  sur¬ 
face  maybe  deformed  spherically  to  account  for  an  incident  wave  of  slight 
spherical  character.  (A  spherical  wave  affects  only  the  coherent  scatter.) 

From  a  measurement  of  scatter  at  normal  incidence  the  surface  power  spectral 
density  and  the  atmospheric  refractive  index  structure  constant  may  be  deter¬ 
mined  provided  the  surface  is  isotropic,  homogeneous,  and  normally  distributed. 

Rough  surface  and  profile  aspects  of  the  theory  are  verified  by  near¬ 
field  measurements  of  scatter  from  titanium  and  aluminum  when  there  is  no 
turbulence  present.  Surface  power  spectral  densities  found  empirically  from 
scatter  at  normal  incidence  agree  with  the  spectra  found  from  proficorder 
tracing. 
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Chapter  I 


INTRODUCTION 

In  the  theoretical  analysis  of  electromagnetic  wave  scattering  from 
rough  surfaces,  there  are  three  factors  in  addition  to  surface  roughness  which 
determine  the  scattering  behavior:  beam  profile,  surface  curvature,  and  prop¬ 
agating  medium  turbulence.  The  rough  surface  statistical  description  is  the 
most  important  aspect  of  the  problem,  and  numerous  previous  investigators 
(most  notably  Beckmann  [1])  have  dwelt  mainly  on  the  scattering  from  zero- 
mean  rough  surfaces  under  spatially  uniform  illumination.  However,  when  tar¬ 
get  dimensions  are  greater  than  the  beam  size,  the  spatial  wave  profile  affects 
the  scattering.  Eckart  [2]  included  a  beam  profile  term,  but  he  restricted  him¬ 
self  to  slightly  rough  surfaces  by  approximating  the  local  normal  on  the  surface 
with  the  normal  direction  of  the  mean  plane.  Also,  surfaces  most  likely  have  a 
mean  structure  that  is  not  a  flat  plane.  Hence  the  general  surface  model  is  a 
composite  sum  of  a  deterministic  and  a  random  component.  Except  when  the 
beam  is  propagated  through  vacuum,  the  medium  will  distort  the  plane  charac¬ 
ter  of  the  wave.  Propagation  in  a  turbulent  medium  is  carefully  examined  by 
Tatarski  [3]  and  Strohbehn  [4]  but  not  in  a  connotation  of  surface  scattering. 

We  assume  the  incident  wave  is  plane,  but  the  field  profile  multiplier  is  a  prod¬ 
uct  of  a  deterministic  term  (the  profile  without  turbulence)  and  a  random  com¬ 
ponent  (which  characterizes  the  medium) . 

The  expression  for  scattered  power  flux  is  derived  from  the  Helmholz 
wave  equation  and  is  averaged  over  ensembles  of  rough  surfaces  and  incident 
waves.  The  final  result  is  a  spatial  Fourier  transform  of  the  product  of  three 
coherences  (or  transfer  function) .  The  first  transfer  function  involves  the  known 
aspects  of  the  problem  —  the  wave  profile  sans  turbulence  and  the  known  sur¬ 
face  component  —  and  is  called  the  DTF  (deterministic  transfer  function) .  The 
second,  the  familiar  optical  (or  modulation)  transfer  function  (OTF^ ,  relates 
the  statistical  properties  of  the  medium.  Finally,  the  rough  surface  transfer 
function  (RSTFl  is  found  from  the  joint  characteristic  function  of  the  surface 
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roughness.  The  advantage  of  the  transfer  function  formulation  is  that  some 
aspects  of  the  problem  may  be  changed  without  affecting  the  others. 

■  !  .  .  I 

The  expression  for  mean  scattered  power  is  applied  to  the  scattering 
of  laser  radiation  from  rough  metallic  surfaces  where  the  medium  is  the  clear 

1  '  II  ' 

atmosphere  (particle  scattering  neglected).  This  application  serves  an  an  ex¬ 
emplary  use  of  the  generalized  result ■,  but  more  important  the  problem  is  the 

,  ■  -  ■  ,  • 

basic  research  in  the  ^tudy  of  images  generated  by  scanning  a  laser  beam 
ovdr  a  target.  In  practice  a  highly-collimated  CO^  or  Nd-YAG  beam  is  scanned 

across  the  target  (in  the  pattern  of  a  Lissajou  or  a  sinusoid,  ,for  example).  If 

*  ,  ( 

the  beam  su.e  is  much  less  than  typical  target  dimensions ,,  an  image  may  be 
created  by  proportionally  displaying  the  received  scattered  radiation  as  in- 

1  (  :  l'  ,  .  , 

tensity  on  a  display  device  suph  as  a  CRT.  Ah  example  of  the  use  of  such  a 

! 

system  is  the  identificatiqn  or  classification  of  aircraft  by  an  air  traffic  con¬ 
troller.  To  study  the  imagery  theoretically,  vye  much  calculate  the  power 

.  i 

scattered  by  the  surface  under  the  beam  at  selected  points  in  time.  To  form 
• '  '  '  ,  , ,  / 

a  suitable  image,  we  expect  surface  roughness  to  give  a  detectable  back- 

I  ’  *  '  ,  | 

scatter  level  even  whqn  the  backscatter  direction  to  the  detector  is  not  near 
'specular.  t  '  ■  >  , 

f  .  .  !  i  ! 

The  theoretical  analysis  was  carried  out  because  of  the  inadequacies 
of  previous  theories  in  predicting  the  propagation  and  scatter  of  a  laser  beam. 

•  :  I  i  f  ! 

Medium  turbulence  (for  example  atmospheric)  becomes  most  significant  when 
the  wavelength  is  short.  Also  beam  profile  and  surface  curvature  ard  signifir 
cant  when  the  beam  is  well-collimated.  Certainly  profile  is  unimportant  when 

I  ! 

a  microwave  radar  beam  is  scattered  from  an  aircraft  but  becomes  significant 
when  the  scattering  is  from  'the  occ  an  surface.  Our  results  apply  at  longer 
wavelengths  such  as  microwave,  UHF,  or  acoustic,  but  some  of  the  four 

.  I  !  I 

aspects  may  be  negligible.  ,  , 

In  the  modeling  of  the  three  system  transfer  functions,  the  laser 

beam  profile  is  the  lowest  order  (TEM  )  circular  mode,  the  deterministic  sur- 

oo 

face  is  quadric,  and  the  log-amplitude  and  phase  fluctuations  of  the  wave  due 
to  turbulence  have  isotropic  Gaussian  statistics.  The  aluminum  and  titanium 


I 


I 


- — 

I 


rough  surface^  studied  are  shown  by  measurement  (tracing  on  a  proficorder  to 
determine  instantaneous  surface  height)  and  by  argument  (using  the  central 
limit  theorem)  to  be  normally  distributed.  From  the  measurement,  the  power 
,  spectral  density  is  modeled  as  a  sum  of  two  Gaussians.  These  two  rough 
surface  components  are  shown  to  be  the  "roughness"  and  "waviness"  compon¬ 
ents  of  the  randomly  rough  surface. 

To  justify  the  use  of  the  models  and  to  verify,  in  part,  the  theory, 
we  present  measurements  of  power  scattered  from  aluminum  and  titanium  air¬ 
plane  surfaces.  Since  the  scattering  is  from  surfaces  with  a  flat  mean  under 
laboratory  conditions  of  no  turbulence,  the  measurements  are  a  verification  of 
the  rough  surface  and  profile  aspects  of  the  theory  and  model.  As  a  by¬ 
product  of  the  theory,  we  learn  how  to  find  the  surface  statistics  from  a 
measurement  of  scattered  power  at  normal  incidence  from  an  isotropic,  homo¬ 
geneous  normally-distributed  surface.  Such  measurements  are  made,  and  the 
surface  height  power  spectral  densities  found  are  compared  with  those  found 
in  proficorder  tracing. 

I 

The  following  sections  include:  a  general  derivation  for  the 
scattered  power,  assuming  the  incident  wave  and  surface  are  known  exactly 
(ch.  II,  sect;  A);  averaging  over  ensembles  of  rough  surfaces  and  turbulent 
waves  (ch.'ll,  sect.  B);  a  suitable  modeling  of  the  transfer  functions  which 
comprise  the  mean  scattered  power  solution  (ch.  Ill,  sect.  A-E);  investiga¬ 
tion  of  special  cases  (ch.  Ill,  sect.  F-G);  and  presentation  of  experimental 
results  and  analysis  (ch.  IV).  Three  related  mathematical  analyses  are  re¬ 
served  for  the  appendix. 
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Chapter  II 

GENERALIZED  THEORETICAL  SCATTERING  ANALYSIS 

This  chapter  presents  the  derivation  of  the  generalized  scattering 
equation  which  includes  the  effects  of  arbitrary  wave  profile,  deterministic 
surface,  rough  surface,  and  the  propagating  medium.  The  result  allows  the 
user  to  specify  his  own  surface  and  wave  properties  for  his  particular  scattering 
problem.  Chapter  III  will  demonstrate  a  suitable  modeling  for  the  laser 
scattering  problem  mentioned  in  the  introduction,  but  for  the  present  all  assump¬ 
tions  are  made  with  the  widest  range  of  application. 


A.  DETERMINISTIC  SCATTERING 

This  section  is  the  derivation  of  the  scattered  power  where  the 
surface  and  incident  wave  are  known  exactly.  The  result  will  be  shown  to  be 
formally  independent  of  incident  polarization.  However,  initially  the  scattered 
electric  field  for  an  incident  field  of  horizontal  polarization  is  derived.  "Hor¬ 
izontal  polarization"  implies  the  incident  E-vector  is  normal  to  the  xz  plane 
of  incidence,  and  the  incident  H-vector  lies  in  the  plane  of  incidence.  (For 
that  plane  and  for  ease  in  understanding  the  subsequent  definitions,  refer  to 
figure  1.)  In  this  section,  Beckmann's  [1]  notation  is  followed  when  possible. 

The  starting  point  of  the  derivation  is  the  scalar  Helmholz  integral 
relation  for  the  scattered  E- field  (E^)  at  a  receiver  point  P: 

E„{P)  =  T-  J7  (  E,  |t-  -  t  -U§>  )  ds  (1) 

2  4  tt  2  3n  fin 

?  E0 

where  in  the  integrand  E„  and  — —  are  the  scattered  electric  field  and  its 

2  o  n  ^  . 

normal  derivative  on  the  scattering  surface  S,  and  i )  and  — -  are  the  Green's 

n 

function  and  its  normal  derivative.  Beckmann  has  the  total  field  at  the  surface 
(E^  +  E^  ,  where  E^  is  the  incident  field)  in  the  integral  rather  than  just  the 
scattered  field.  However,  Stakgold  [5]  shows  the  two  expressions  are  equal 
by  proving 


J7  <E 

S 


1  ?n 


*  )  ds  5  0. 

a  n 


(2) 
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Figure  1.  Surface  geometry.  (Only  one  surface  point,  Ps ,  is  depicted  to 
reduce  complexity.) 

The  appropriate  Green's  function  is 

*  =  (3) 

where  k  is  the  wave  number  and  R'  is  the  variable  distance  from  P  to  a  surface 

point  P  .  If  the  receiver  point  P  is  in  the  far-field,  only  the  first  two  terms 
s 

of  the  power  series  for  R'  need  be  kept: 

kR'  «  kR  -  k„  •  r  (4) 

o  2 

where  Rq  is  the  (constant)  range,  k^  is  the  propagation  vector  directed  from 
the  origin  to  P,  and 

r  =  x  +  y  +  C  (*,y)  z  .  (5) 


(The  hat  (  A  )  above  the  letter  indicates  a  unit  vector:  for  example  z  =  z  /  |  z  I .) 


The  familiar  £  is  chosen  as  the  surface  height  function.  From  (3)  and  (4)  we 


find 


(6) 


exp  (ikRn)  .  — 

ill  «  — ^ exp  (-ik2  *  r)  . 

o 


In  the  far-  field  limit,  the  vector  directed  from  P  to  P  is  assumed  to  have  the 
_  s 

unit  vector  k  /  k  . 

4 

The  incident  electric  field  is  assumed  harmonic  and  nearly  plane: 


Ej  (r  ,t)  =  U  (r)  exp  (  ik  .  r  -  iiot  )  .  (8) 

(In  keeping  with  Beckmann  notation,  sub-1  implies  an  incident  quantity; 
sub-2,  scattered.  As  usual,  the  temporal  frequency  is  tn,  and  t  is  time.) 

The  wave  profile  factor  U  (r)  is  complex  in  general  and  includes  (a)  the 
deviation  of  the  E-field  from  a  pure  plane  wave  (due  to  scintillation),  and 
(b)  the  original  spatial  profile  at  the  wave  source.  In  equation  (8),  k^  is 
a  propagation  vector  directed  from  the  far-field  source  point  Q  to  the  origin 
of  the  scattering  system. 

The  tangent  plane  approximation  is  made  as  a  preliminary  to  a 
Kirchhoff  boundary  condition.  The  surface  must  have  large  local  radii  of 
curvature  with  respect  to  X,  the  wavelength,  for  the  approximation  to  be  rea¬ 
sonable.  Moreover,  the  slope  along  the  surface  must  be  small  (less  than  a 
few  degrees)  so  that  the  incident  polarization  on  the  plane  as  a  whole  is 
equivalent  to  the  polarization  on  the  individual  tangent  planes.  Both  Leader 
r  6]  and  Valenzuela  r7]  have  shown  that  depolarization  is  an  effect  of 
multiple  scattering;  but  since  small  slope  and  large  radius  of  curvature  are 
assumed  here,  the  surface  roughness  is  slight  to  moderate,  and  multiple 
scattering  will  be  insignificant.  Finally,  since  the  simple  plane-wave 
reflection  problem  is  solved  to  establish  the  boundary  conditions  on  a  single 
tangent  plane,  we  are  justified  in  using  the  Fresnel  horizontal  reflection 
coefficient  R  in  expressing  the  scattered  E-field  at  the  surface  as  R  times 
the  incident  field  there: 


(9) 


What  to  use  for  the  Fresnel  coefficient  is  a  subject  of  continuing  argument; 
however,  for  present  purposes,  we  assume  R  is  not  a  function  of  the  integration 
variables  of  (1).  This  is  at  least  valid  for  infinitely  conductive  surfaces;  and 
as  an  approximation  for  other  surfaces,  an  average  R  taken  over  the  surface 
for  a  given  incident  angle  might  be  sufficient.  Error  for  this  type  of  average 
increases  with  increasing  roughness  of  the  surface,  but  once  again  the  tangent- 
plane  approximation  restricts  us  to  moderately  rough  surfaces. 

The  normal  derivative  boundary  condition  is  also  found  from  the 
tangent-plane  reflection  problem: 


E1R 


(-ik, 


.  J  1  MK 

*  +  l; 


do) 


In  most  circumstances  the  second  term  in  the  parentheses  in  (10)  is  negligible. 
Mathematically  we  require 


1 


k  U  cos  0 


1 


*U 

az 


«  1 


(11) 


and  that  the  surface  slopes  ( £  ,  f  )  be  much  less  than  unity.  This 

x  y 

says  grazing  incidence  is  avoided,  and  the  gradient  of  U  normal  to  the  mean 
surface  is  smaller  than  kU.  There  is  an  example  of  a  verification  of  (11) 
in  appendix  I  wherein  the  inequality  is  checked  for  the  models  presented  in 
chapter  III.  Satisfying  (11)  also  implies  there  is  only  a  slight  variation  in 
U  in  the  z-direction  over  distances  the  order  of  the  wavelength  and  hence  of 
the  composite  surface.  If  the  Maclaurin  series  in  z  for  U  (r)  is  written: 


U  (r)  =  U  (x,y,0)  +  ? 

5  z 


z  + 


(11a) 


z=0 


then  by  inequality  (11)  only  the  first  term  in  the  series  need  be  kept: 


U  (r)  «  U  (x , y) .  (12) 

After  putting  (6),  (7),  (8),  (9),  (10),  and  (12)  into  (1),  dropping  the 
time-harmonic  factor,  and  performing  vector- algebra  manipulation: 
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E,(P)  -  S4  eXBP(‘kRo)  IT  Ueiv'r[-C  dxdy  (13) 

2  4  ttR  x  x  y  y  z 


where  v=k,  -k„=vx+vy  +  vz  (14) 

1  2  x  y  z 

=  k  [  (  sin  0  -sin  0  cos  0  )x  -  sin  P  sin  0  y-(cos  0  +cos  0  )z  ] . 

*  Ct  O  Ck  J  1  4 

(Subscripts  x  and  y,  when  appearing  on  £  or  U,  imply  those  partial  derivatives.) 

Equation  (13)  is  the  first  reasonable  form  for  use  when  U  and  £  are  exactly 

known,  but  without  further  approximation  the  integral  is  revised  via  partial 
integration  such  that  the  derivatives  appear  on  U  rather  than  £: 

CD  -  - 

E  (P)  =  C  R  [  f  [U - (v  U  +  v  U  )  ]  eiV  r  dx  dy  (15) 

2  1  J  J  u  v  x  x  y  y  J 


where 


ik  exp(ikR0)  ,  vx  ,  ~  \ 

—  ,VP  (cosV^  Sin  V 


2  —  —  2 

and  v  =  v*v  =  2k  (l-  sin  0  sin  R  cos  0  +  cos  0  cos  P  ) .  (17 

1  “  O  1  4 

Reduction  of  (15)  by  eliminating  the  second  term  in  the  bracket  requires  (for 
each  (x,y)  and  (©  ,  ©2  .  G3)): 


1E1  - 

2k  (|Ux|  +  |Uy|) 


>>  1. 


loosely,  (18)  is  satisfied  if  (a)  there  are  no  edges  or  sharp  boundaries 
(shadowing  is  neglected);  (b)  k  is  large  (this  is  consistent  with  the  local 
radius  of  curvature  being  much  larger  than  X);  and  (c)  if  0=0,  then 
^1=  (grazing  incidence  is  avoided  if  the  observer  at  P  "looks  at" 

the  specularly  scattered  component).  If  (18)  is  indeed  satisfied,  then  the 
expression  for  scattered  E-field  of  horizontal  polarization  for  known  U  and  £ 


00  iv  £  i(v  x  +  v  y) 


E2(v)  =  C  R"  JJ  U  e  z  e 


dx  dy 


which  is  written  to  emphasize  the  Fourier  kernel.  Like  inequality  (11),  the 
validity  of  inequality  (18)  is  checked  in  appendix  I  for  the  specific  field 
profile  and  medium  turbulence  model  of  chapter  III. 
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The  above  result  is  also  applicable  to  scattering  of  a  vertically 
polarized  incident  electric  field: 

=  U(r)  [  x  cos  @1  +  2  sin  ]  exp( ik^  •  r  -iu>t)  (20) 

=  E,  x  +  E,  z 
xx  z 


or  an  H-field  of  either  polarization.  The  new  electric  field  boundary  conditions 


Ej  =  R  (-E.  x+E,  ±)\  (21) 

^|S  x  z  |g 

a"d  4r  TJ>  <El  s‘El  l)  <22> 

9n  1  Uon1x  ZS 

+  ^  4- 

where  R  is  the  vertical  Fresnel  coefficient.  (For  the  forms  for  R  and  R  , 
see  Beckmann  Tl]  /  p.  21.)  Because  of  the  similarities  between  (9)  and  (21), 
and  (10)  and  (22),  the  steps  leading  to  the  scattering  equation  (19)  are  found 
to  be  applicable  to  each  of  the  two  scalar  components  of  the  scattered  field 
(termed  E2  and  E2  )  where  the  incident  field  is  vertically  polarized.  The 

X  Z  - 

(scalar)  Euclidean  norm  (  jEyT  +  E22  )  of  the  scattered  field  is  the  desired 

X  z 

quantity  since  it  is  the  field  observed  by  a  receiver  at  P  looking  along  ko- 
That  scalar  scattered  field  is  also  given  by  (19)  if  we  employ  inequality  (11) 
again  and  let  R  -*R+.  The  scattered  magnetic  fields  are  then  found  by  duality. 
For  example,  to  derive  the  scattered  H-field  due  to  an  incident  H  lying  normal 
to  the  incident  plene,  replace  the  E  in  the  derivation  of  scattered  E-field  for 
horizontal  polarization  by  H,  and  let  R  -R+. 

The  real  scattered  power  (flux)  is  now  easily  written  for  either 
vertical  or  horizontal  polarization: 


P  =  |  En  H 
2  2  2  2 


_  i  v  •  A  _ 


=  ~Y  I  ci  1  R  JJ  Gd(vz;  A)  e  1  d 


where  R  is  the  reflection  coefficient  of  the  appropriate  polarization,  and  Yq 
is  the  admittance  of  free  space  .  The  function  G  in  (23)  is  an  unnormalized 
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mutual  coherence: 


is  a  spatial  distance,  and 

d  A  =  d  Ax  d  Ay.  (25) 

In  a  coherence  in  the  "bra-ket"  form  in  (24),  the  convention  is  to  take  (x,y) 

as  the  subscript  1,  and  (x  +  &x,  y  +  Ay)  as  sub-2.  Not  only  is  the  functional 

form  of  G  ,  a  coherence,  but  the  Fourier  transform  in  (23)  of  G,  is  a  power 
d  d 

spectral  density  function.  Furthermore,  note  that  in  general  a  rovariance 
B  ( A)  which  is  defined  by 

B  (A)  =  <  (VJ-<Vl>)*  (V2  -<V2»  >  (25a) 

is  in  turn  related  to  the  coherence  by  rewriting  (25a): 

B(A)  =  <V*V2>  -<V*><V2>  ,  (25b) 

and  that  a  correlation  C  ( A )  is  defined  by 

C(A)  =  -Hi!  .  (25c) 

B  (0 ) 


In  general,  the  incident  polarization  is  neither  purely  vertical  nor 
horizontal  but  is  arbitrary.  When  superposition  is  used  on  the  incident  E 
and  H  fields  to  represent  an  arbitrary  polarization,  one  finds  that  the  relation 
(23)  is  still  applicable.  Suppose  is  rotated  through  an  angle  a  from  the 

incident  plane  about  the  propagation  vector  k  (as  depicted  in  figure  2).  Then 

2  * 
the  coefficient  R  needs  only  to  be  redefined: 


(26) 


2  +  2  2 
R  =  (R  cos  a)  +  (R  sin  a)  . 

However,  the  integrals  (23)  and  (24)  are  unaltered.  The  preceding  polar¬ 
ization  analysis  indeed  shows  that  the  term  "scalar  scattering"  is  still  ap¬ 
propriate  even  when  the  incident  fields  are  of  arbitrary  polarization.  In 

fact  the  question  of  incident  polarization  can  be  ignored  since  it  is  obviously 

2 

a  matter  of  choosing  the  proper  form  for  R  .  From  now  on,  the  case  of  a  perfect 
2 

reflector  (R  =  1)  is  taken  without  significant  loss  of  generality. 


z 


Figure  2 .  Rotation  from  a  vertical  incident  polarization  (dashed  lines) 
to  an  arbitrary  one.  Note  each  vector  E1  and  can  now  be 
split  into  vertically  and  horizontally  polarized  components. 
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B.  ENSEMBLE  AVERAGING 


Since  very  often  in  scattering  the  surface  is  not  exactly  known,  the 

scattered  power  as  averaged  over  an  ensemble  of  random  surfaces  is  needed. 

The  actual  surface  model  chosen  is  a  sum  of  a  randomly  rough  (C^)  and  a 

deterministic  surface  (£  ,): 

d 

c=cr  +  ed.  (27) 

The  rough  surface  is  described  by  its  second  order  statistics  and  is  assumed 

to  be  statistically  homogeneous  (stationary)  on  the  area  under  the  beam.  The 

joint  probability  density  of  the  random  variable  £  is  W  (z  ,z  ;  Ax,  Ay),  and 

r  l  m 

the  joint  characteristic  function  is  the  double  Fourier  transform: 


x(v1<v2;  A)  =  JJw(z1#z2;a)  e  ^  V*  Z*  +  ^  dz]L  dz2 

—  CO  ^ 

=  <  e  i(Vl  Zl  +  VsZa)>  . 


(28) 


(The  reader  should  be  careful  not  to  confuse  the  statistical  characteristic 
function  (28)  with  the  deterministic  coherence  function  (24)  even  though  the 
bracket  notation  is  the  same.)  The  surface  ensemble  average  is  now  taken 
on  P2  given  by  (23): 

<p2>  =  \  ICJ2  JJGd(vz;A)x(vz;A)eiV-L'A  dA  (29) 


where  is  as  in  (24)  except  now  £  -*  £  ,  and 

x(vz;  A)  =  x<V  “vz  ;  • 


(30) 


Homogeneity  of  allowed  the  y  to  be  extracted  from  the  integral  defining 

G  .  If  the  surface  model  chosen  were  a  more  general  higher  order  composite: 
d  n 


C  56  C  + 


(31) 


i=l 

and  if  y  is  the  joint  characteristic  function  of  £  ,  then  in  (29): 
1  i 

ns 

y  ( v  ;  A )  =  H  y,  (v  A  ) 
z  i=l  *i  z 


(32) 
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if  {  Cr  }  is  a  mutually  independent  set. 

Next,  rather  than  assuming  exact  knowledge  of  the  incoming  wav/e, 
an  average  over  an  ensemble  of  incoming  waves  U  is  taken.  The  field  profile 
is  rewritten  as  the  product: 

U  =  C  z  (33) 


where  P  is  the  deterministic  character  of  the  incident  field,  and  Z  is  a 
homogeneous  random  variable  whose  statistics  are  assumed  known.  The 
random  variable  Z  is  the  output  of  a  turbulent  medium  where  the  input  is  a 
spatially  uniform  illumination  of  unity  amplitude.  (Hence  Z  is  itself 
unitless.)  A  second  averaging  is  then  performed  on  the  mean  scattered  power 
over  the  ensemble  of  incident  waves: 

<<p2>>  =  !  CjJ2  JjGd(v2;A)  x(vz;A)  Cq(A)  e  ^  A  dA  (34) 

-a? 

<n> 

where  C  ( A )  =  <  Z,  Z„  )  is  a  statistical  mutual  coherence  (unnormalized),  and 
o  1  Z 

the  G(j  is  still  defined  by  (24)  except  now  U-  p  and  Q  -  Q  .  In  obtaining 
(34),  the  homogeneous  property  of  the  medium  (and  hence  of  Z)  is  used  to 
extract  Cq  from  the  integral  defining  G^.  The  model  for  the  wave  profile  (33) 
is  especially  chosen  to  yield  this  decoupling,  since  U  itself  is  non-homo- 
geneous. 

Relation  (34)  is  the  most  general  form  for  mean  scattered  power  flux 

to  be  obtained  here.  It  sets  forth  the  four  contributors  to  the  problem  —  the 

incident  wave,  transmitting  medium  turbulence,  the  deterministic  surface, 

and  the  rough  surface.  It  quite  obviously  reduces  to  the  earlier,  less  general 

scattering  equations  (23)  and  (29).  This  extends  Beckmann's  result  in  that  his 

involves  only  the  rough  surface.  Three  special  cases  are  notable:  (a)  if  the 

rough  surface  goes  to  zero,  x  -  1;  (b)  if  there  is  no  turbulence,  C  -* 1;  and 

o 

(c)  if  the  mean  surface  is  a  flat  plane  (  C  =  0),  G  ,  becomes  the  mutual 

d  d 

coherence  which  Eckart  [2]  found  as  part  of  his  result: 
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K^=o 


=  [J  P (x,y)  e  (x  +  Ax,  y  +  Ay)  dxicty 


=  <  *1  V 


The  above  integral  may  also  be  termed  the  incident  wave  coherence. 

Interpretation  of  (34)  is  best  done  via  a  systems  transfer  function 
viewpoint.  Since  (34)  is  a  Fourier  transform  (where  v  is  spatial  frequency), 

i 

each  component  of  the  integrand  (except  the  kernel)  can  be  interpreted  as  a 


transfer  function,  as  exemplified  by  the  "system"  of  figure  3.  In  usual  termi- 

i  , 

nology  C  is  the  optical  (or  modulation)  transfer  function  (OTF) .  Let  G  ,  and 
o  d 

X  be  defined  as  the  deterministic  and  rough  surface  transfer  function^  (DTF  and 


RSTF),  respectively.  The  "input"  to  the  system  is  the  constant  factor 

Y  !C  I  .  The  amplitude  of  the  incident  wave  is  in  the  DTF  G_.  Because 
o .  1 '  d 

of  the  decoupled  nature  of  the  result,  certain  aspects  of  the  problem  (such 
as  the  rough  surface)  can  be  revised  without  varying  others  (deterministic 
Components  and  the  turbulence) . 


TRANSMITTING 
MEDIUM  RESPONSE 


Gd  (vz;  A) 
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>« 
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KNOWN  SURFACE 
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ROUGH  SURFACE 
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Figure  3.  System  transfer  function  model  for  the  mean  scattered  power  as 
given  by  equation  (34). 
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Chapter  III 

MODELING  OF  TRANSFER  FUNCTIONS 

The  general  solution  (34)  is  now  applied  to  the  scattering  of  laser 
radiation  from  rough  curved  surfaces  where  the  propagating  medium  is  the  clear 
atmosphere .  This  serves  as  an  exemplary  use  of  the  result  and  more  important 
gives  the  desired  theoretical  information  on  scattering  from  the  titanium  and 
aluminum  airplane  surfaces.  Initially  the  rough  surface  statistics  are  inves¬ 
tigated  to  find  the  RSTF;  then  the  OTF  and  DTF  are  found.  Although  these 
models  are  optimal  for  the  problem  at  hand,  they  are  also  applicable  in  numer¬ 
ous  parallel  situations.  The  advantage  of  the  formulation  (34)  is  that  inap¬ 
propriate  transfer  functions  may  be  remodeled  without  disturbing  the  others. 

For  example,  suppose  after  calculating  the  scatter  from  the  titanium  surface 
in  the  atmosphere,  we  desire  the  scatter  from  the  same  surface  under  water. 
Then  tie  new  OTF  would  be  determined  but  the  remainder  of  the  analysis 

i  «  • 

would  remain  intact. 

I  ' 

A.  ROUGH  SURFACE  TRANSFER  FUNCTION 

The  aluminum  studied  was  taken  from  a  scrapped  USAF  F-84  tail 
section.  Because  the  part  had  been  in  extensive  use,  it  is  more  typical  of 
aircraft  in  use  today  than  an  unused  sample  would  be.  Aluminum  is  generally 
used  in  sub-Mach  2.5  aircraft,  and  titanium  in  Mach  2.5  and  up.  The  titanium 
studied  is  a  sample  used  in  stress  testing.  Since  titanium  does  not  corrode  as 
aluminum  does,  it  too  is  typical  of  present-day  aircraft  materials.  Other  sam¬ 
ples  which  were  not  investigated  were  (a)  a  painted  one,  since  presently  paint 
of  military  aircraft  is  classified  and  hence  unavailable;  (b)  composite  materials 
such  as  Boron/Epoxy  and  Graphite/Epoxy  which  are  formed  by  a  process  which 
does  not  obey  the  central  limit  theorem.  In  fact,  the  composites  exhibit  a 
nearly  periodic  surface.  Moreover,  neither  composite  nor  painted  surfaces 
have  a  near-unity  reflection  coefficient. 


To  establish  the  RSTF,  the  second-order  statistics  of  the  rough 
surface  must  be  found  by  directly  probing  the  surface.  Since  the  surfaces 
under  study  are  optically  rough,  the  proper  instrument  is  a  proficorder.  A 
proficorder  (Bendix  trade-name)  drags  a  stylus  slowly  across  a  surface  and 
produces  an  output  proportional  to  surface  height.  The  analog  output  is 
recorded  and  converted  to  digital  form  for  computer  analysis.  In  our  case,  the 
recording  was  made  on  a  strip-chart  and  the  conversion  was  manual.  Traces 
were  made  on  the  surfaces  at  various  locations  to  check  for  homogeneity 
(stationarity);  moreover  at  some  locations,  traces  at  90°  and  45°  to  each 
other  were  made  to  check  for  isotropy.  A  total  of  six  traces  with  number  of 
points  ranging  from  128  to  512  were  taken  for  the  aluminum  to  give  a  total 
number  of  1920  data  points  for  analysis.  There  were  four  titanium  traces, 
each  with  128  data  points.  The  sample  period  was  0.2  mil;  hence  trace 
lengths  ranged  from  25.6  to  102.4  mil  in  length.  The  first  step  in  computer 
analysis  is  the  removal  of  a  linear  trend  and  mean  to  minimize  surface  orient¬ 
ation  errors.  The  output  of  the  program  includes  probability  density  (p.d.f.), 
cumulative  distribution  (c.d.f.),  and  power  spectral  density  (p.s.d.f.). 

1.  Probability  Density. 

Williamson  [8],  [9]  has  shown  that  to  remarkable  accuracy  most 
randomly-prepared  metallic  surfaces  possess  a  Gaussian  probability  density. 
For  surfaces  where  the  preparation  technique  is  a  cumulative  result  of  a  large 
number  of  repeated  events  occurring  randomly  over  the  surface,  the  central 
limit  theorem  applies  and  insures  a  nearly  normal  probability  density  [8]. 
Such  preparation  must  completely  define  the  surface  roughness;  it  cannot 
merely  reshape  maxima.  Hence  sanding,  sand  blasting,  corrosion,  and 
rolling  of  sheet  metal  may  generate  a  normal  density;  whereas  light  polish¬ 
ing  and  buffing  may  not.  Moreover,  if  the  effect  of  these  processes  is 
distributed  uniformly  over  the  surfaces,  the  statistics  will  be  homogeneous. 
And  if  the  process  operates  without  "preferred  direction",  the  statistics  are 
isotropic.  Hence  the  above-mentioned  processes  should  generate  isotropic, 
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homogeneous  normal  statistics.  (An  example  of  a  non-isotropic,  but  homo¬ 
geneous  random  surface  is  the  ocean.  Here,  the  wind  direction  is  the  pre¬ 
ferred  direction.)  Rather  than  relying  completely  on  this  heuristic  approach, 
the  surfaces  are  checked  for  p.d.f. ,  isotropy,  and  homogeneity. 

The  probability  density  of  the  aluminum  and  titanium  surfaces 
was  obtained  by  the  histogram  method  suggested  by  Jenkins  and  Watts  [10], 

To  see  how  close  to  a  Gaussian  the  p.d.f.  is,  the  c.d.f.  was  found  by  a 
cumulative  summation  of  the  p.d.f.  and  was  plotted  on  probability  paper. 

The  composite  c.d.f.  for  titanium  in  figure  4  and  for  aluminum  in  figure  5 
(indicated  by  "total")  closely  follow  a  straight  line;  so  these  surfaces  are  nearly 
normally  distributed  .  That  straight  line  is  the  c.d.f.  corresponding  to  a 
Gaussian  p.d.f.  which  has  the  same  variance  and  mean  as  found  for  the 
sample.  In  figure  4,  the  abscissa  is  normalized  to  the  standard  deviation; 
however  the  abscissa  in  5  is  not  normalized.  The  two  c.d.f.'s  in  figure  4 
labeled  "roughness"  and  "waviness"  are  explained  in  part  2  of  this  section. 

Not  only  were  the  composite  p.d.f.'s  nearly  Gaussian,  but  the  c.d.f.  for 
each  surface  trace  (not  shown)  closely  followed  the  straight  line  as  well. 

The  individual  standard  deviations  are  given  in  the  captions  of  figures  6 
and  7.  The  properly-weighted  average  standard  deviations  are  :  aluminum 
0.067  p,  ;  titanium  0.095^.  From  this  information  we  can  conclude  the 
assumption  of  stationarity  to  first  order  is  reasonable.  We  still  need  to  show 
second  order  homogeneity  is  valid. 

2 .  Power  Spectral  Density. 

Since  the  surfaces  can  be  suitably  represented  by  a  normal  p.d.f., 

all  that  remains  in  completely  defining  the  surface  second-order  statistics  is 

a  measurement  of  the  power  spectral  density,  P  .  The  p.s.d.f.  is  obtained 

s 

through  the  following  steps:  (a)  taking  the  (fast)  Fourier  discrete  series 
expansion  of  the  data,  (b)  applying  a  Hanning  window  as  a  three-point  con¬ 
volver,  (c)  taking  the  modulus-square  of  the  resultant  sequence  as  a  raw 
estimate  of  power,  and  (d)  averaging  over  bands  of  frequency  to  improve  con¬ 
fidence  limits.  (The  final  limits,  or  error  bars,  were  evaluated  from  the  table 
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provided  by  Blackman  and  Tukey  [ll].)  The  individual  p.s.d.f.'s  for 

each  of  the  four  titanium  traces  are  shown  in  figure  6;  and  figure  7  has  the 

six  aluminum  p.s.d.f.'s.  The  titanium  spectra  display  the  expected  isotropic 

and  homogeneous  behavior;  however  the  aluminum  is  at  best  only  isotropic 

and  locally  homogeneous  in  that  traces  labeled  1,2,3  are  at  one  location  on 

the  surface  and  traces  4,5,6  are  one  inch  away.  The  poor  behavior  of  the 

aluminum  is  due  to  non-uniform  corrosion  of  the  surface.  Hence 

for  the  aluminum,  an  averaged  power  spectrum  is  used  to  represent  the  entire 

surface.  Since  the  titanium  appears  stationary  (to  second  order),  averaging 

of  titanium  P  improves  confidence  limits.  (Titanium  averaged  p.s.d.f. 
s 

is  in  figure  8;  aluminum,  figure  9.)  When  scattering  measurements  are  made 
and  compared  with  theoretical  curves  which  use  these  averaged  spectra, 
we  will  observe  a  good  match  which  confirms  the  necessity  of  taking  the 
average  for  aluminum. 

We  notice  that  there  are  two  regions  in  which  the  power  is  located. 
The  low  frequency  power  is  termed  the  "waviness"  component  because  of  the 
long  wavelengths.  Waviness  is  a  result  of  the  rolling  process  that  formed  the 
metal.  Similarly,  the  wide-band  component  is  termed  "roughness" .  The 
origin  of  the  roughness  component  is  partly  the  process  which  formed  the  metal 
originally  and  partly  the  processes  that  act  on  the  surface  after  formation  (e.g. 
conrosion  in  the  case  of  aluminum) .  Because  of  the  multiple-application  of 
events  to  create  the  random  surface  structure  in  surface  formation  and  in  the 
post- format  ion  processes,  we  expect  the  roughness  and  waviness  components 
each  to  be  normally  distributed.  Each  surface  trace  was  then  filtered  into  its 
two  components,  and  each  was  investigated  for  c.d.f.  Once  again  the  exper¬ 
imental  results  confirm  the  expected  Gaussian  p.d.f.  The  c.d.f. ' s  of  the 
individual  titanium  components  are  shown  in  figure  4  wherein  the  abscissa  of 
each  is  normalized  to  the  individual  standard  deviation. 

Since  the  two  surface  components  appear  to  be  created  by  inde¬ 
pendent  processes  (even  in  sheet  rolling  the  small-scale  fluctuations  are  not 
related  to  large-scale),  the  components  are  themselves  independent.  Hence 
the  sum  of  the  individual  power  spectra  equals  the  total  spectrum,  and 
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POWER  SPECTR 
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Figure  6. Individual  titanium  surface  height  power  spectral  density  versus 
wavenumber  v  .  (Each  trace  has  128  data  points.  Individual 
standard  deviations  (in  p,)  are  0.099  (#1),  0.102  (#2),  0.084 
(#3),  and  0.094  (#4).  Nyquist  frequency  is  0.0984  p,-*.) 
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Figure  7 .Individual  aluminum  surface  height  power  spectral  density  versus 
wavenumber  v.  Confidence  limits  are  not  shown  due  to  non- 
stationarity.  (Individual  standard  deviations  (in  p)  followed  by 
number  of  data  points  are  0,072,  512  (#1);  0.062,  128  (#2);  0.125, 
256  (#3);  0.063,  256  (#4);  0.049,  256  (#5);  0.077,  512  (#6). 
Nyquist  frequency  is  0.0984  p  \) 
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the  surface  characteristic  function  can  be  written  as  the  product  of  individ¬ 
ual  characteristic  functions.  For  the  individual  power  spectra  P  ,  a  Gaussian 

si 

model  is  closely  fitted  to  the  data  (see  figures  8  and  9): 

P.  (v)  =  Ps  (0)  exp  [-(v/v0  )2]  (36) 

fai  i  i 

where  v  is  spatial  frequency  and  vQ  is  the  e-folding  frequency  of  the 
Gaussian.  Consequently,  the  individual  surface  autocorrelation  functions 
are  the  Fourier-cosine  transforms  of  P_  (v)  : 


C_  (T)  =  f  cos(2ttvt  )  Pq  (v)  dv  .  (36a) 

Si  -»  i 

Hence  Cg  (t)  =  exp  [-(t/t^)2]  (37) 

where  =  (  Ax)2  +  (Ay)2,  (38) 

_  ]  2  — 

t  =  (ttv  )  ,  and  o0  =  P„  (0)  v  Jr*  is  the  individual  surface  variance. 

°i  °i  si  si  °i 

3  .  Form  of  the  RSTF. 

In  general,  the  RSTF  is  the  joint  characteristic  function  of  the 

surface  roughness  x(v  ,v  ;A)  evaluated  along  the  diagonal  v=v  ;  v  =-v  . 

1  L  1  z  z  z 

But  since  the  composite  surface  given  by  (31)  is  desired,  where  Qr ^  is  the 

"waviness"  component,  and  Cr  is  the  "roughness"  component,  the  RSTF 

2 

(by  (32)  )  is  the  product  of  the  respective  characteristic  functions  evaluated 
along  the  diagonal.  For  each  roughness  component  Qr  the  RSTF  is 

T)  =  e^P  [-a2  v 2  ( 1  -  Cs  ( t )  ) ]  .  (39) 

\  i  Z  1 

Extension  to  higher  order  surface  composites  (  ng>  2)  is  straightforward. 
Depending  on  the  size  of  the  normalized  surface  variance 

gl  =  vz  °si  (40) 

with  respect  to  unity,  the  characteristic  function  x*  may  be  written  two 
ways: 

Case  1:  For  g^  ^  1 ,  the  surface  is  slightly  to  moderately  rough,  and  xA  in 
(39)  is  expanded  in  an  exponential  series  which  in  general  appears  as: 


(41) 


*  m 

2  v  2 

6  ~  L  m!  * 

m=0 

The  plot  of  figure  10  gives  the  absolute  error  ( e)  in  approximating  the 
infinite  series  (41)  by  the  series  which  is  truncated  at  m  ,  as  a  function  of 

*  _5 

z.  For  example  when  z  =  0.1,  m^=  3  is  sufficient  for  e£  10  .  So 

ml  m  2 

X  (v  ;  t)  «  e  9iV  e  m^To^  .  (42) 

i  z  Li  mj 

m=0 

Case  2:  For  g^  »1,  the  surface  is  very  rough,  and  the  characteristic  function 

Xi  is  predominantly  affected  by  Cs^  near  unity  (or  t  near  zero) .  Hence  only 

the  first  two  terms  of  the  Maclaurin  series  for  Cq  are  retained  and 

si 

Xi (vz;t)  » exp  [-gj  (t/t^)2]  .  (43) 

This  case  should  be  used  with  caution,  however,  since  for  very  rough  sur¬ 
faces  the  tangent  plane  approximation  weakens. 

4  .  Checking  Boundary  Conditions. 

Finally,  the  initial  assumptions  of  small  slope  and  large  radius  of 
curvature  are  checked.  The  mean  slopes  ( J 2  as  / t0 from  eq.  (106)) 
are  found  not  to  exceed  0.3  degrees  for  the  waviness  component  and  1.0 
degrees  for  the  roughness  component  .  Hence  the  small  slope  assumption 
is  valid.  To  check  for  large  radius  of  curvature,  the  one-dimensional  cur¬ 
vature  (  K )  is  found  point-wise  on  the  computer  for  the  samples  using: 

K  (x)  *  [ 1  +  (&■)  ]  •  (43a) 

(Curvature  is  the  inverse  of  the  radius  of  curvature . )  The  probability 
density  of  K  is  found  empirically,  and  the  standard  deviation  is  found 
from  that  density.  Moments  of  the  density  of  K  are  more  meaningful  than 
those  of  l/K  since  in  general  the  latter  has  a  bimodal  density  and  the  former 
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a  nearly  monotonic  density.  The  small  curvature  assumption  is  said  to  be 

i  -3  -l 

valid  if  \«  1  /a^r.  In  particular,  Oj,  ^  6.7  x  10  p  ;  hence  \«  150p. 
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Figure  10.  Summation  limit  versus  argument  z  such  that 

i  mi 
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£  e  .  For  determining  the  number 


,  of  terms  needed  to  represent  the  diffuse  component  of  scatter 

*  :  I  , 

when  gi  «  1 . 
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B„  OPTICAL  TRANSFER  FUNCTION 

In  this  section  an  OTF  for  the  clear  atmosphere  is  derived  in  a  form 

I 

suitable  for  the  integration  of  (34) .  Such  a  solution  is  not  only  applicable  to 
the  scattering  problem  but  more  generally  supplies  the  loss  in  coherence  of  a 
plane  wave  in  the  optical  to  infrared  wavelengths  propagated  through  the 
atmosphere.  An  alternate  form  for  the  OTF  derived  by  Hufnagel  and  Stanley[12l 

M  (p )  =  exp  f-c  p5^3]  (44) 


(where  c  is  some  constant) , is  perhaps  accurate  but  prohibits  closed-form 
integration  of(34)  since  (44)  is  not  analytic  at  the  origin.  Rather,  an 
approximate  Gaussian  model: 

CQ(p)  =  exp  [  -c'p2]  (45) 

(c'  a  constant)  is  derived  herein  from  the  initial  steps  provided  by  Berann.31, 
T 1 4 1  and  Strohbehh  T  4  ],  T 1 5 1. 

It  is  inappropriate  here  to  reiterate  Tatarski1  s  work  [3]  on  wave 
propagation  in  a  turbulent  medium.  For  our  purposes,  the  most  useful  quan¬ 
tity  presented  by  Tatarski  is  the  three-dimensional  spectral  density  of  the 
index  of  refraction  fluctuations  as  a  function  of  wavenumber: 


Vk) 


-11/3 


H  <  H  <  X 

o  m 


assuming  local  isotropy  and  homogeneity.  In  (46),  h  =  5.92 /l  and  h  =  2tt/L 

m  o  o  o 

where  l  and  Lq  are  the  inner  and  outer  scales  of  turbulence,  respectively. 

The  outer  scale  (L  )  is  usually  taken  [41  as  the  constant  height  above  ground, 
and  lQ~  1  to  10  mm.  Strohbehn  [4]  has  improved  the  form  of  by  rewriting 

the  function  to  remove  the  singularity  at  the  origin: 

1  1  /  '  2  2  "11/6 

or  [  1  +  x  L  ]  kSk 

o  m 

<M*)  =  (47) 

J 


where 


C2  L11/3 

n  o 


5.41  tt" 


(  k  L  »  1  ) 
m  o 


;  t 

and  C  is  the  refractive  index  structure  constant.  Typical  values  for  C/ 

u  r 


i « -.TtmT*  W«»'  ivrw»»yi 


ir-'.a ti  >  lyi. ■  ’<*  '>0 


(for  L  <  a  few  hundred  meters)  given  by  Brown  [16]  range  from  6  x  10“^ 

-2/3  -15 

m  '  for  the  strongest  daytime  turbulence  (L  <  a  few  meters)  to  6  x  10 

-2/3  0 

m  for  very  weak  turbulence  at  dawn  or  dusk.  Actually,  Strohbehn  used 

2 

a  Gaussian  taper  function  expC-U/n  )  ]  to  turn  off  $  at  k  =  k  .  However, 

m  n  m 

due  to  the  low  value  of  $  at  h  ,  form  (47)  is  suitable  here.  As  the  wave- 

n  m 

number  h  decreases  for  r<  kq,  the  spectral  density  loses  validity  for  two 
reasons:  (a)  The  model  presented  by  Strohbehn  is  only  supposed  to  reduce 
the  mathematical  problems  by  removing  the  singularity  at  the  origin.  It  does 
not  necessarily  fit  experimental  data  for  small  k.  (b)  For  such  small  h, 
the  spectrum  actually  loses  meaning  since  one  may  not  assume  homogeneity 
over  very  large  separations  in  the  corresponding  structure  function  [3]. 

The  OTF  is  found  in  terms  of  by  assuming  the  statistics  of 
the  wave  amplitude  (  |z|  )  are  log  normal  and  the  wave  phase  (arg(Z)  )  are 
normal.  Experimental  measurement  [  17] , [18]  and  theoretical  analysis  [3], 
[4]  using  the  central  limit  theorem  have  validated  those  assumptions.  Beran 
[13], [14]  and  Strohbehn  [4], [15]  have  derived  the  OTF  in  terms  of  the  wave 
structure  function  (D)  by  somewhat  independent  procedures: 

C  (A?  ,  ATI)  =  exp[-D(AF  ,  AT|)/2]  (48) 

o 

where  (AF  ,  ATi)  are  distances  in  the  "incident  system"  (F  ,  T))  which  is 

normal  to  the  incident  propagation  vector  k^.  (The  wave  structure  function 

is  the  sum  of  Tatarski's  log-amplitude  and  phase  structure  functions.)  If, 

in  addition  to  statistical  homogeneity,  isotropy  is  assumed,  then  C  =  C  (p) 

o  o 


where 


p2=  (A?)2  +  (ATI)2. 


The  wave  structure  function  is  related  in  general  to  the  wave  covariance 
(B  )  by 

D(  p )  =  2  (BJO)  -  Bu(p)  ).  ( 

If  the  wave  correlation  function  C^(  p)  =  Bu(p)/  B^O),  where  6^(0)  is  the 
wave  variance,  then 

C  (p)  =  exp[-B  (0)  (1  -  C  (p)  )  ].  (! 

o  u  u 
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Note  the  formal  similarity  between  (51)  and  the  RSTF  (39).  In  (51),  B  (0)  is 

2  2  2  U 

a  wave  variance  normalized  to  X  ,  and  in  (39)  as  vz  is  a  normalized  surface 
variance.  The  reason  for  the  similarity  is  that  a  joint  characteristic  function 
involving  normal  random  variables  is  obtained  as  the  RSTF;  whereas  a  coher¬ 
ence  of  log-normal  variables  is  the  OTF. 

The  expression  for  the  OTF  in  (43)  and  (51)  is  for  an  infinite  plane 
wave  Z.  The  results  of  this  section  are  applicable  to  the  scattering  problem 
only  if  the  separation  U  =£Z  in  (33)  is  valid.  In  an  independent  approach  based 
on  Whitman  and  Beran  [20],  the  decoupling  is  not  assumed,  and  precisely  the 
same  scattering  result  (34)  is  derived.  The  reader  can  verify  this  by  using  an 
arbitrary  boundary  condition  rather  than  Whitman's  (3)  and  then  evaluating  the 
resulting  coherence  corresponding  to  his  (22)  at  frequency  p,  =  0. 

The  correlation  is  written  as  a  Fourier- Bessel  transform  of  the 
spectrum  $  (h): 

n  2 

Cu(  p)  =  J  Iq(hp)  4>n(n)  k  dn  (52) 


where  Jq  is  the  Bessel  function  of  first  kind,  zero  order.  At  this  point  we 
note  that  since  for  small  k  loses  validity,  CM  p)  for  large  p  (>Lq)  is 
likewise  questionable.  Through  straightforward  integration  the  wave  var¬ 
iance  is  found  to  be: 

00 

B  (0)  =  4  tt2  k2  R  f  $  (  h  )  h  d K 
u  o  J  n 

o 

=  0.786  k2  R  C2  L5/3.  (53) 

o  n  o 

Equation  (52)  must  be  integrated  numerically,  except  for  small  p«t  when 

o 

c(p)  «.  i  -  4ra  h  )1/3  -43(  p/L„>2-  <54> 

u  4  o  m  b  o 


(Equations  (52),  (53),  and  (54)  all  assume  h  L  »1.)  Unfortunately,  the 

mo  2 

absolute  error  in  using  (54)  is  less  than  or  equal  to  (  pn  )  per  cent.  Since 

m 

the  beam  size  is  much  greater  than  l  ,  form  (54)  is  essentially  us’eless.  So 

o 

computer  integration  is  used  to  give  a  close  approximation  to  Cu(  p): 
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Cu(0)  «  exp  [-1.125  (Bu(0)  )1/4  (p/LQ)2  1,  3 . 3(B°(0)  )5/8  (55) 

where  B°(0)  =  B  (0)  L  is  not  a  function  of  k  and  only  depends  on  L 
u  u  o  m  o 

to  the  extent  that  does.  The  steps  in  deriving  (55)  are  outlined  in 

appendix  II.  For  h  <3.3  (B°(0)  ,  equation  (54)  is  more  appropriate 

m  u 

than  (55);  however,  (54)  will  be  rarely  needed  except  at  long  ranges  or  very 

short  wavelengths.  Furthermore,  for  B  (0)  »  1  (easily  achived),  C  (  p)  in 

u  o 

(51)  is  affected  predominately  by  CM  p)  near  unity  (i.e. ,  for  p  near  zero). 
Hence,  only  the  first  two  terms  of  the  Maclaurin  series  for  C^(  p)  are  kept 
in  writing  the  OTF  in  the  desired  Gaussian  form: 

C  (  p)  »  exp  [  -(  p/p  )2]  (56) 

o  o 

where  p  =  0.945  (B°(0)  f5/8  L  _l/24.  (57) 

o  u  o 


For  L  less  than  a  thousand  meters,  the  e-folding  length  p  can  be  approx- 
o  o 

imated  as 

p  »  1.1  (  k2  R  C2)  _5/8,  (58) 

o  on 

-1/24 

which  is  not  quite  dimensionally  correct  since  the  Lq  factor  has  been 
dropped.  Since  C  (  p)  for  p  near  zero  is  all  that  is  needed  in  writing 
the  OTF,  the  questionable  character  of  CMp)  for  large  p  is  of  no  concern. 
Equation  (58)  says  that  the  OTF  of  the  atmosphere  is  only  a  function  of 
wavelength,  range,  and  the  atmospheric  structure  constant. 

It  is  now  interesting  to  compare  the  above  OTF  with  Hufnagel  and 
Stanley' s  OTF  [  1 2  ]  : 

M(  p )  =  exp  [ -1.455  k2  Rq  C2  p5/3  ]  .  (59) 

A  graphical  comparison  of  (56)  and  (59)  reveals  only  slight  discrepancies, 
and  indeed  (56)  is  more  suited  for  the  closed-form  integration  of  the  mean 
power  given  by  (34).  In  that  integration,  the  independent  variable  p  will 
be  written  in  terms  of  the  scattering  coordinate  system  (x,y)  rather  than  the 
incident  system  (f  ,T|)  using  the  isotropic  distance  p  where 

p2  =  (  Ax  cos  e^2  +  (Ay)2.  (60) 
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Another  OTF  which  may  be  of  value  to  the  reader  was  derived  by 
Yura  [19]  for  transmission  through  ocean  water.  He  shows  two  mechanisms 

_3 

give  rise  to  small-angle  (less  than  10  rad)  forward  light  scattering: 
suspended  biological  particles  with  permittivity  approximately  that  of  water, 
and  large-scale  (with  respect  to  X)  refractive  index  fluctuations.  The  OTF 
for  the  latter  mechanism 

M(p)  =  exp  r-(p/pQ)2]  (60a) 

where  p  =  2.76  (  k2  R  C2  L1/6  (60b) 

o  o  n  o 

is  of  course  quite  comparable  to  the  atmospheric  OTF  in  (56)  and  (57).  Yura's 


notation  has  been  changed  for  comparison  with  the  present  author’s  results. 
Note  that  Yura's  variance  of  index  of  refraction  fluctuations  is  An2  = 


0.526  L2/3  C2 
o  n 


Initially  it  appears  that  the  e-folding  length  of  the  OTF 


is  larger  in  water  than  in  the  air  —  a  result  which  is  not  expected.  However, 

2  -  8  -  2  /3 

this  is  not  the  case,  since  in  water  is  typically  about  3  x  10  m  , 

much  larger  than  in  air.  Notice  it  is  a  trivial  matter  to  change  from  a  cal¬ 
culation  of  scattering  from  a  target  in  the  a  nosphere  to  scattering  from  the 

same  target  underwater.  All  that  changes  is  the  value  of  the  "constant"  p  . 

o 

Because  of  the  formal  similarity  of  the  OTF  in  (51)  with  the  RSTF 
in  (39),  a  rough  surface  with  appropriate  statistical  characteristics  can  be 
used  to  simulate  the  turbulent  atmosphere.  Rather  than  transmit  the  wave 
through  the  atmosphere,  the  beam  would  be  reflected  at  normal  incidence  from 
a  rough  surface  (or  transmitted  through  a  random-phase  transparency).  The 
normally-distributed  surface  would  have  height  variance  B^(0)  /(4  k  ),  and 
the  e-folding  length  of  the  Gaussian  autocorrelation  would  be  pQ.  The  sur¬ 
face  would  typically  have  a  much  larger  variance  and  correlation  length  than 
the  aluminum  and  titanium  surfaces  under  study;  hence  the  roughness  could  be 
formed  by  peening.  Moreover,  the  mean  slope  /Bu(0)  / k  pq)  (from  (106)) 

would  not  exceed  about  1.0°  (using  X=  10.6  p,  ,  R  =  L  =  1000.  m,  and 

2  -14  -2/3  0  ° 

C  =  6  x  10  m  ),  and  so  the  tangent-plane  approximation  would  still 
n 

be  valid.  This  is  of  course  expected  because  the  propagating  medium  cannot 
distort  the  phase  of  the  wave  so  as  to  violate  the  original  boundary  conditions. 
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C.  DETERMINISTIC  TRANSFER  FUNCTION  .  „ 

iv  Cd 

The  DTF  G  .  is  given  by  (24)  where  V  =  £  e  .  For  the  wave1 
d 

profile  £,  the  lowest  order  (TEM  )  circular  output  mode  of  the  laser  (assum- 

oo 

ing  uniform  spatial  coupling  at  the  output  mirror  )  is  taken  in  the  scattering 
coordinate  system: 

e(x,y)  =  £q  exp  {-f(x  cos  e^2  +  y2  ]/l2  )  (61) 

where  l  is  the  e-folding  length  of  the  electric  (magnetic)  field  profile  and 

£  is  the  real  amplitude  factor.  The  e-folding  length  (£')  of  the  incident 
o 

power  is  related  to  l  by  l  =  1‘  JY  . 

The  deterministic  surface  model  is  a  quadric.  In  a  coordinate 
system  (u,v)  which  is  rotated  through  an  arbitrary  angle  7)  from  (x,y)  but 
in  the  same  plane  as  (x,  y): 

£d(u'v)  =  '  (  W  +TT  )  •  (62) 

U  V 

The  transformation 

(u,v)  =  (x  cos  71  -  y  sin  Tj ,  y  cos  7)  +  x  sin  7]  )  (63) 

is  needed  to  write  {^(x,y).  Relation  (62)  is  the  most  general  needed  for 
a  second-order  surface  since  terms  of  lower  order  can  be  removed  by  rotation 
about  the  x  or  y  axes.  A  rotation  about  the  x-axis  can  be  done  arbitrarily 
since  that  merely  redefines  the  incident  polarization  which  has  no  effect 
here  since  R  =  ±  1 .  That  rotation  would  not  be  proper  in  the  case  of  a  non¬ 
perfect  conductor,  and  hence  first-order  terms  would  be  needed.  A  rotation 
about  the  y-axis  can  always  be  freely  done  since  the  incident  angle  6^  is 
with  respect  to  the  normal  to  the  surface  mean  at  the  origin  of  the  scattering 
system.  However,  rotation  about  the  z-axis  is  not  allowed  since  the  orienta¬ 
tion  of  this  known  surface  (set  by  T] )  with  respect  to  the  incident  (xz)  plane 
cannot  be  changed  at  will. 

If  the  exact  surface  can  be  adequately  represented  by  terms  to 
second  order  in  its  Maclaurin  series,  then  the  quadric  surface  (62)  is 
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L 


suitable.  Such  surfaces  include  a  small  area  on  a  sphere,  an  ellipsoid, 
and  a  circular  cylinder.  In  fact  p  (  or  p  )  is  the  radius  of  curvature  of  £ 


u  'v  ad 

taken  along  the  v-axis  (or  u-axis)  and  then  evaluated  at  the  origin  u  =  0 


(or  v  =  0).  Moreover,  the  form  in  (62)  will  allow  for  closed-form  integration 
for  the  DTF  and  the  final  power.  Higher  order  terms  of  the  expansion  force 
computer  integration. 

Integration  for  the  DTF  yields  the  (anticipated)  Gaussian  form: 

[-■f?  ~  T~  (ei(«)2+2S2ixiy+e3(iy)Z)] 


G  ,(v  ;A)  =  - 2- 

d  z  cos  6 


exp 


1 


(64) 


2  2  2 

where  A  -  tt  l  ft  /2  is  a  pure  constant  formed  for  conservation  of  power, 
o  o 

Namely  £  Y  A2  is  the  total  power  in  the  incident  wave  neglecting  turbu- 
o  o 


lence.  Also 


Pj  =  sec201  ( 


2 

cos  y 


2  2 
sin  1)  ■ 


u 


V 


)  ♦  sin2*  cos2!,  (±-±) 


u 


V 


P =  -  Sin  y  cos  y  (  — - -  ) 

2  K  pu  pv  ’ 


r  2 

'  sec  6 


,( 


2  .  2 
cos  7)  +  sin  n 


u 


)  +  ( 


.  2 
Sin  Vi 


cos  y 


u 


/  _1 


and  p.  =  sec20  (  — - —  ^  sin2,n  cos2Ti  +  ( 

j  1  v  p  p  /  \ 


2 

sin  V 


2  2 
+  °°LA)  . 

p  / 


(65) 


Special  cases  including  (a)  a  parabolic  cylinder  (p^-  ®  ),  (b)  a  paraboloid  of 


revolution  (o  =  p  ),  and  (c)  a  flat  plate  ( p  -*  »  ,  p  -*  ®  )  may  be  easily  worked 
u  v  u  v 


out  if  desired.  For  example  if  the  mean  surface  is  a  flat  plane,  the  DTF 
becomes  the  incident  wave  coherence: 


A2 

o 


Cd=0 


cos  0 


exp  (-fp)  • 


(65a) 


Hence  the  e-folding  length  of  the  incident  wave  coherence  is  l,  of  the 
incident  field  profile  is  l,  and  of  the  incident  power  flux  is  1/  JT  .  The 
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radii  of  curvature  p^  and  p^  may  be  either  positive  or  negative;  however, 
to  justify  the  use  of  the  tangent-plane  approximation,  |puJ  and  |pv|  must 
be  much  greater  than  the  wavelength  and  the  beam  size.  Hence  the  scattering 
result  (34)  may  not  give  the  correct  mean  scattered  power  for  every  portion  of 
the  target  surface  —  especially  if  the  beam  size  is  large  (e.g.,  one-fifth  the 
length  of  the  target),  and  the  surface  curvature  is  large  (e.g., on  a  wing- tip). 

D.  FINAL  INTEGRATED  RESULT 

The  general  result  (34)  with  the  three  transfer  functions  (OTF 
(56,57),RSTF  (32,39),  and  DTF  (64)  )  substituted  in  is  integrated  to  yield 
the  expression  for  mean  scattered  power  in  closed  form.  Since  Case  1  above 
(g.  **  1;  i=l  ,2)  is  better  suited  to  the  original  tangent  plane  approximation, 
and  since  the  surfaces  being  studied  are  slightly  to  moderately  rough,  only 
the  Case  1  expansion  is  necessary  for  each  component  of  the  rough  surface. 
After  integration  of  (34): 


A2  |C.  |2  Y  .  .  "J1  ni  m  n 

_  o  i  o  “(g,  +ga)v  v  g-.  ^ 


« P2» 


2  cos  9. 


m!  n! 


m=0  n=0 


G(Y;vx,vy)  (66) 


where  G  is  defined  as  a  generalized  two-dimensional  Caussian  integral: 

00 

G(y;u,v)  =  JJ  exp  [i(u  l  +  vll)  -  (  y1^2+y2?T1  +  YgT)2 )  ]  d?dn 


4  YlVY| 


-u2  y3+  uvy2-v»Yi' 

4  Y  Y  -  Ya 
Y1  3  Y2 


(y1>0'  4yiy3>Y2K 


In  particular,  the  arguments  y  for  (66)  are: 

/  1  .  1  \  2  .  .  m  ,  n 

Yl  =  (  •?  +TP)  cos  VlT  +  + 


2  ,2  Q 
v  l  B, 
z  1 
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i  •  mi.*  n. , 


Y2  = 


2  #  2  O 

Vz  1  82 


and 


V3  = 


=  (_!,_  +  1 

V  p3  2-t,3  /  T  3 


n 


2  . 2q 

V  B3 
8 


(68) 


One  further  solution  is  written  for  the  convenience  of  the  reader  —  that  of 

a  single  component  very  rough  surface  (Case  2,  g  »1  ;  n  =  1): 

s 


«V>  = 


Aa 

o 


icj' 


2  cos  6 


1 


G{Y;  vx'V 


(69) 


where  y  is  as  in  (68),  except  let  n  -*  0,  m-*g,  and  tq  -  tq  . 

Solutions  of  the  form  (66)  and  (69)  are  useful  in  computer  program¬ 
ming  the  solution  to  the  scattering  problem  but  give  little  insight  as  to  the 
mechanisms  affecting  the  scattering  behavior.  Next,  to  reduce  the  complexity 
of  the  solution  (66),  typical  values  of  the  parameters  of  y  in  (68)  are  inves¬ 
tigated. 


E.  NUMERICAL  RESULTS 

The  rough  surface  parameters  were  determined  in  the  analysis  out¬ 
lined  in  section  A  of  this  chapter.  For  all  the  surfaces  examined,  0.04  ^ 
Og^  0.1  p,  (i  =  1,2);  hence  the  roughness  parameter  gi  ^  2  for  X  ?>  1.0  p.  . 
This  justifies  the  use  of  only  Case  1  in  the  previous  section,  since  the 
Nd-YAG  (X  =  1.06  p.)  and  the  CO^  (X  =  10.6  p)  lasers  are  the  most  promising 
imaging  sources.  Except  for  very  long  ranges,  very  short  wavelengths,  or 
very  small  illuminated  areas; 


,  i =  1,2. 


(In  fact,  the  definition  of  waviness  and  roughness  components  implies 

t  3  »  t3  ,  but  that  assumption  is  not  needed  here.)  Notice  that  if  the  wave 
°1  °2 

were  propagated  underwater,  p  would  be  much  smaller  than  in  the  atmosphere, 

o 
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and  the  above  assumption  would  not  be  applicable.  (Refer  back  to  (60b)  and 

(58).)  Furthermore,  a  flat  mean  (C  =  0  =  8.;  i=l,2,3)  is  taken  to  reduce  the 

d  i 

complexity  and  to  increase  our  understanding  of  the  result. 

The  solution  (66)  is  then  separable  into  four  significant  terms: 


;>’>>  -  $■  |C‘„|3Y°  e-k'+9*){AexPr--^2£il((  nv2)] 

2  2  cos  9  l  L  4  tt  \\  cos  9,  /  y  /  J 
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+  To  "  I  WiV  exp  f-  Jox  \  /  <4m>  ]  +  to22"  L  Vi hr  exp[  ‘V  v1Z/(4n)] 


m=l 
mi  n: 


ni 

+Tvr 

n=l 


-  exPr.  v2/(4(  a+  A  ))  1} 

LL  ml  +  JL)  P|  1  y  Jn! 

m=l  n=l  V  tT  t  2  J  1  2 

°1  °2  (70) 

2  -  - 

where  v  =  v  •  v  ,  and  an  effective  area  of  illumination  is  defined  as 
xxi 


A  =  "[««81(IiT  +  4r)] 


-1 


(71) 


In  the  special  case  of  no  turbulence  (  p  -  <*>)  and  a  flat  mean  (  C,=  0),  the 

O  Q 

effective  area  is  the  total  integral  of  the  incident  coherence  (A2  2  {®tt/cos29, ) 

o  1 

divided  by  the  peak  incident  coherence  (A2  /  cos  9,).  (The  incident  coherence 

o  1 

is  given  in  general  by  (35)  and  for  a  Gaussian  profile  by  (65a).)  Alterna¬ 
tively  the  area  A  is  the  area  of  an  ellipse  with  semi-axes  equal  to  the  e-folding 
lengths  in  the  Ax  and  Ay  directions  of  the  incident  coherence  ( ./?  1/  cos  0^ 
and  J1 1 ,  respectively). 

The  first  term  in  (70)  is  defined  as  the  coherent  (or  specular) 
scatter  since  it  is  the  observed  scatter  from  a  perfectly  flat  plate  (i.e. , 
gx  =  g2  =  0).  The  second  term  is  the  primary  diffuse  scatter  due  to  the  surface 
waviness  component;  third  is  the  primary  diffuse  term  due  to  the  surface  rough¬ 
ness  component;  and  the  double  summation  is  a  remaining  contribution  due  to 
all  surface  roughness.  If  g  «  0.25  (where  g  =  gi  +  ga  ),  then  the  double 
summation  in  (70)  is  negligible.  In  that  case  we  can  let  nx  =  mx  =  1  since 


_3 

the  absolute  truncation  error  will  not  exceed  10  (obtained  from  figure  10 
for  g  =  0.2).  Hence,  if  g  «  0.25  : 


«p2» 


2  cos  0, 


11  To,  9i  exp['T0i^  /4]  +  "V,  g2  exp["ro2v12  / 4  ]}  •  <72> 


Equation  (72)  is  appropriate  for  the  CO  laser  wavelength,  10.6 p,  ,  since 

Ct 

g  0.02  for  the  titanium  and  aluminum  under  examination.  However,  if 
the  Nd-YAG  laser  is  used,  X  =  1.06p,  g^2.,  and  all  the  terms  in  (70)  must 

be  retained.  (The  summations  could  terminate  at  n  =m  =8  for  an  error  less 

-3  11 

than  10  . ) 

Note  that  in(71)  an  increase  in  turbulence  (p  decreasing)  implies 

o 

a  decrease  in  effective  area  A.  Consequently  the  coherent  component  widens , 

and  its  peak  decreases.  Most  important,  the  coherent  term  is  affected  by 

both  the  wave  profile  (via  l)  and  the  turbulence  (via  p  ) ;  whereas  the  diffuse 

o 

terms  depend  only  on  surface  parameters  (  t_  and  ac  ).  If  the  turbulence 

°i  i 

had  been  much  stronger  (say  for  a  wave  propagated  underwater),  then  pQ 
would  have  influenced  the  diffuse  terms  rather  than  the  coherent.  The  expres¬ 
sion  for  such  a  case  will  not  be  shown  since  it  is  taken  care  of  by  (66).  When 
the  surface  becomes  rougher  (by  increasing  g),  power  is  coupled  from  the 
coherent  component  to  the  diffuse.  That  decrease  in  the  coherent  component 
is  accounted  for  by  the  exp  f-g  ]  factor  which  multiplies  the  entire  expression. 
It  is  observed  that  variable  (  and  hence  a  variable  deterministic  surface  C^) 
will  affect  either  the  diffuse  or  the  coherent  terms,  or  both,  for  typical  0,.  If 
such  values  are  lo  be  included,  then  one  should  either  use  the  general  expres¬ 
sion  (66)  or,  based  on  the  actual  values  of  p. ,  decide  whether  the  known 
surface  affects  the  diffuse  or  the  coherent  terms  in  (70). 

To  give  further  feel  for  the  scattering  expression  (70),  sample 
scattering  patterns  at  the  CO  laser  wavelength  are  given  in  figure  11  for 

id 

titanium  and  figure  12  for  aluminum.  The  condition  of  turbulence  is  weak 
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Figure  11.  Polar  plots  of  mean  scattered  power  (in  db)  for  0  =0,  0 

3 


constant,  6  varied;  normalized  to  power  at  0  =0  =0  =0 
2  12  3 


Parameters  for 'titanium  found  in  surface  tracings 


Figure  12.  Polar  plots  of  mean  scattered  power.  Same  conditions  as 
figure  11  but  for  aluminum  parameters. 


2  -15  -2/3 

(C  =  6  x  10  m  ,  R  =  250  m)  so  that  the  individual  components  of 
i  n  o 

scatter  may  be  discerned.  The  polar  plots  of  scattered  power  are  in  decibels 
where  the  reference  power  is  the  backscattered  power  at  normal  incidence. 
(Hence  the  norm  is  not  a  function  of  (0  ,0  ,0  )  but  does  depend  on  the  sys- 

1  b  J 

tern  parameters  l  ,  p  ,t0  ,  and  os  .)  For  each  incident  angle  (8  =  0°, 30°  , 60°) , 
the  widest  and  next-to-widest  lobes  are  the  diffuse  scatter  due  to  roughness 
and  waviness  respectively.  The  extremely  narrow  lobe  is  the  coherent  com¬ 
ponent.  In  practice  wider  coherent  components  are  encountered  because 
of  stronger  turbulence  and  because  of  an  incident  spherical  wave.  All 
reasonable  degrees  of  turbulence  are  taken  care  of  by  (70);  how  to  account 
for  an  incoming  spherical  wave  is  the  subject  of  the  next  section,  part  1. 

<  ■  I 

F.  SPECIAL  CASES 

In' this  sectiop  we  give  further  interpretation  of  the  scattering 
result  (70)  by  investigating  special  cases.  Actually  the  first  case  extends 
(70)  by  introducing  incident  waves  with  slight  spherical  character.  Second, 
in  the  special  case  of  zero  mean  {£d=0),  equation  (70)  reduces  to  Beckmann’ s 
result.  Third,  the  direct  backscatter  case  is  analyzed  to  show  more  closely 
the  way  roughness  parameters  affect  scattering  behavior.  Finally,  scatter 
from  a  perfectly  smooth  surface  is  presented.  All  cases  are  tied  closely  toge¬ 
ther,  and  the  reader  will  note  some  conclusions  apply  to  many  special  cases. 

1 .  Incident  Spherical  Waves. 

As  pointed  out  in  the  examples  of  scatter  in  figures  11  and 
12,  the  coherent  component  of  scatter  is  often  much  wider  than  predicted  for 
an  incident  plane-wave.  The  reason  in  part  is  that  the  incident  wave  must 
always  possess  some  spherical  (or  second-order)  curvature.  To  account  for 
this,  we  do  not  have  to  go  back  to  the  Helmholz  integral  and  rewrite  the 
expression  for  incident  fields.  Rather,  the  deterministic  surface  is  deformed 
spherically  with  radius  of  curvature  twice  that  of  the  radius  of  the  incident 
wave-fronts.  If  the  incident  divergence  angle  is  a  (between  e-folding  points 
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'T^jrrum 


rr.i-iiw'WT'PFifT-W- 


of  the  incident  field  ),  then  Tl=  0,  o  =  4  1/ a  .  and  0=0/  cos  0  .  Hence 
2~  v  _o  u  v  1 

Pj=  62=  Co/ (4 <•)]  ,  02=  0  is  substituted  iny,  eq.  (68).  As  a  result,  the 

revised  coherent  term  (first  in  the  brace  of  (70)  )  becomes 


=  «r-s-(44)] 


yvjvj 

^  l 

4 

where 
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?)  +  T 
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z  o 
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(73) 


provided 


2  /  ,  /  »2 
°o  <<  (  XAo  } 


i=  1,2. 


To  demonstrate  the  significance  of  including  an  incident  spherical  wave, 

suppose  p  =10  cm,  /  =  1  cm,  a  =1  mrad,  and  \  =  10  5  m.  Then  1/(2  f.3)= 
3  _2  0  2  ^4_2 

5x10  m  and  (v  a)  /128  s;  1 . 25x  10  m  .In  this  typical  example,  turbulence 
z  o 

is  negligible,  but  the  spherical  wave  is  as  important  as  the  wave  profile 
contribution. 


2 .  Normalized  Power,  Compared  With  Beckmann. 

The  normalized  expression  for  mean  scattered  power  is  shown  to 

agree  with  Beckmann's  result  [  1  ]  in  the  special  case  of  a  zero-mean  rough 

surface  (C  =  0),  a  single-component  rough  surface  model  (n  =  1),  and  a 
d  s 

perfectly  plane  incident  wave  (a  =  0).  Normalization  is  with  respect  to  back- 

o 

scatter  at  normal  incidence  from  a  perfectly  smooth  surface.  The  area  A 

defined  by  (71)  is  especially  constructed  to  correspond  to  Beckmann's  area 

4  L  L  ,  where  2  L  and  2  L  are  the  sides  of  the  uniformly  illuminated  rec- 
x  y  x  y 

tangular  surface.  Since  turbulence  is  lumped  into  the  effective  area  A,  it 
is  not  necessary  to  require  no  turbulence  in  the  comparison.  The  normalized 
expression  is: 
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„2  -gf  r  A  cos  9,  ,2  I 

P  =  F3  e  {  exp  L  -  TJ— 1  V  J 


m 


m! 

-£nz  e*pt-T#vV(4mn} 

m=l 


(74) 


where 


,2  _  f  Vx  N*  .  2 

V  \  cos  0j_  /  Vy 


and 


g  =  (v,  <0  * 

z  s 

V 


f3  =  [ 1  -  t:  tan  »i  ]  • 


Beckmann's  expression  (eq.  (57),  p.  88  in  f  1  ]  )  can  be  achieved  if  the  notation 
is  revised  by 

To  -  T  • 

2  2 

V  -  V  , 

i  xy 

-  «  , 

_  2  r  A  cos  9,  ,  2  l  2 

and  F3  exP  L  4^  1  v  J  -  P0  • 

The  last  expression  is  the  normalized  coherent  component  of  scatter  (the 
scatter  from  a  perfectly  smooth  surface  replacing  the  rough  one) . 

3 .  Direct  Backscatter . 

When  an  image  is  formed  by  scanning  the  source  beam  over  the  tar¬ 
get,  a  large  amount  of  backscattered  power  is  needed  for  all  incident  angles. 

To  investigate  how  the  power  in  direct  backscatter  can  be  increased,  the 
power  as  given  by  (70)  (at  0  =  -0  ;  0  =0  J  is  normalized  to  the  backscattered 

£»  1  J 

power  at  normal  incidence  (a  pure  number)  by: 


^P2^B 


(6,-0,O) 


/  «P2» 


(0,0,0) 


(75) 


To  simplify  the  analysis,  we  take  only  a  single -surface  model  (n  =1).  Then 

s 
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(<P0))D  =  (sec40)  exp  T(  2  k  a  sin  0  )2  1 

Co  S 


r  r  Ak2  -  _s  _  .T1  (a  2kc°s  0)2m 

{  eXp[-  —  sin  e  tan  9  j+  L  - ^ -  exp. 


(T0k  sin  0)' 
m 


m=l 


]}• 

(76) 


If  kTQ  is  decreased,  the  surface  roughens  because  the  peaks  come  closer 

together  relative  to  the  wavelength.  Hence  the  diffuse  scatter  widens 

Since  we  have  no  control  over  t  ,  we  must  increase  the  wavelength  to  Increase 

o 

backscattered  power  at  the  large  angles,  0.  Unfortunately,  a  decrease  in 

kT  causes  a  decrease  in  the  peak  (at  0=0)  and  total  power  of  the  diffuse 
o 

component.  So  there  is  a  trade-off  between  desired  width  and  peak.  In 

2 

addition,  an  increase  in  surface  variance  a  will  increase  diffuse  back- 

s 

scatter  overall.  An  analysis  similar  to  the  above  can  be  performed  if  the 
imagery  is  effected  at  rone  bistatic  angle  (i.e.  receiver  not  located  at  laser 
sourer!  yet  monostat.c  Ti>agery  gives  us  the  required  information  on  the  effect 
of  surface  parameters  «he  -'iffuse  scatter.  Sample  backscatter  plots  using 
the  aluminum  and  titarium  parameters  are  provided  in  figure  13.  These  plots 
and  those  o.  figures  11  and  12  are  generated  on  a  digital  computer  using  the 
composite  result  (70).  because  of  the  low  level  power  in  the  diffuse  scatter, 
it  is  readily  obvious  that  receiver  electronics,  laser  power,  and  detector 
sensitivity  will  have  to  be  quite  good  to  detect  backscattered  power  for  angles 
significantly  different  from  the  specular. 


4 .  Scatter  From  a_ Perfectly  Smooth  Surface. 

In  order  to  provide  a  heuristic  verification  of  the  coherent  com¬ 
ponent  in  (70),  the  scatter  from  a  perfectly  smooth  surface  for  vertical  in¬ 
cidence,  when  there  is  no  turbulence  present,  is  investigated.  The  normalized 

coherent  component  of  (74)  for  p  — •  <»  is : 

o 

p  *  exp  [-(  0  /0  )2  ]  (77) 

c  2  20 

where  V)*  (78) 
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Figure  13.  Polar  plots  of  mean  backscattered  power  (in  db)  for  6  =  o, 

J 

6  =  =  -02  varied;  normalized  to  power  at  =  9„  =  0^  =  0. 


is  the  e-folding  angle  (  0?  «tt/2  assumed  ). 

0 

profile: 


Suppose  divergence  dominates 


1 

TV 


« 


k3  a* 
32 


(79) 


Then  the  e-folding  angle  is  the  same  as  the  e-folding  angle  of  the  incident 
power: 

""  (80) 


e2  =  S- 
2o  JQ 


(Recall  a  is  the  total  divergence  angle  measured  between  e-folding  points 
on  the  incident  field  profile.)  Now  suppose  profile  dominates  divergence: 


k5  at 
32 


(81) 

(82) 


This  compares  favorably  with  the  half-angle  divergence  1 . 36/(kje)  of  a  uniform 
beam  diffracted  by  a  round  pin-hole  of  diameter  4 Jl  l.  Hence  for  both  the 
special  cases  of  the  coherent  scatter,  the  solution  reduces  to  well-known 
results  of  diffraction  theory. 


G.  THE  INVERSE  PROBLEM 
1 .  Scatter  at  Normal  Incidence . 

The  subsequent  analysis  demonstrates  how  to  determine  the  surface 
power  spectral  density  and  the  atmospheric  refractive  index  structure  constant 
from  a  measurement  of  scattered  power  where  the  wave  is  normally  incident  on 
the  surface  (0^=0).  All  the  assumptions  used  in  the  modeling  apply  here 
except  we  assume  no  knowledge  of  the  surface  power  spectrum.  One  addi¬ 
tional  restriction  is  that  the  surface  is  only  slightly  rough: 

k2  a2  «  1  ,  (83) 

s 

2 

where  is  the  rough  surface  variance.  Recall  that  of  all  the  aspects  of 
modeling,  a  form  for  P  was  never  theoretically  derived.  Rather,  a  sum  of 


two  Gauss ians  was  determined  experimentally  to  be  a  good  approximation. 
Once  again,  the  assumptions  of  statistical  isotropy  and  homogeneity  for  a 
normally  distributed  surface  are  applied.  If  those  assumptions  are  valid, 
then  a  measurement  of  the  surface  power  spectrum  gives  a  complete  second- 
order  statistical  description  of  the  surface.  Since  this  method  determines 
properties  of  the  surface  and  the  medium  from  scatter  measurements,  the 
analysis  is  referred  to  as  the  "inverse  problem . "  Not  only  are  these  results 
valuable  in  studying  surfaces  at  large  ranges,  but  also  in  getting  a  quick 
statistical  description  of  a  surface  in  a  laboratory.  There  is  less  error  in¬ 
volved  in  this  passive  probe  since  most  surface  tracing  devices  physically 
deform  the  surface. 

The  revised  expression  for  normalized  mean  scattered  power 

(for  0=0  )  is  once  again  a  sum  of  a  coherent  (p  )  and  a  diffuse  term  (p 
1  c  d 


P(02>  =  Pc  +  pd 


=  exp 


A'  v2 
4 


v  3  2  n 
zo 

A' 


[*  B  (  t  )  J  (vt)  t  dt  (84) 
0S  ° 


where 


L'  “  "  [  27* 


kfa? 

32 


v  =  k  sin0. 


and 


v2  =  k2  (  1  +  cos  0  )2 
zo  2 


The  surface  covariance  B  (i )  =  a2  C  (t)  is  assumed  to  drop  to  a  value  much 

s  s  s _ 

less  than  a2  in  a  distance  t  «  ,/A'  /tt  ,  and  that  the  e-folding  angle  of 
s 

the  coherent  component  is  much  less  than  tt/2.  The  Fourier- Bessel  transform 
in  (84)  is  achieved  because  of  vertical  incidence  and  because  of  the  isotropy 
of  the  rough  surface  and  the  atmosphere.  Normal  incidence  makes  the  three 
transfer  functions  dependent  on  t  rather  than  some  being  functions  of  p. 

The  coherent  term  in  (84)  gives  the  information  on  the  atmosphere 
by  relating  pQ,  the  e-folding  length  of  the  OTF,  to  the  measurable  coherent 
component  of  scatter.  The  e-folding  angle  of  the  coherent  component  ©2 
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is  related  to  A'  through 


Rearranging  equation  (85)  we  find: 


Finally,  the  refractive  index  structure  constant  C  is  related  to  p  via 

n  o 

equation  (58)  or 


<p„/  i.i  r8/s 


(87) 


Hence  the  strength  of  turbulence  can  be  measured  by  measuring  the  e-foldir.g 

angle  of  the  coherent  scatter  (assuming  we  know  R  ,  k,  a  ,  and  t ). 

o  o 

From  the  diffuse  component  of  scatter  p  in  (84),  we  can  find  the 

d 

surface  power  spectral  density  (and  hence  the  variance  and  autocorrelation). 


Since  the  surface  two-dimensional  power  spectral  density  Vg  is  the  Fourier- 
Bessel  transform  of  the  covariance  (see  appendix  III ),  then  by  direct  com¬ 
parison  with  diffuse  component  in(84): 

A'  p.(v) 

Vv)-  TtTv5  ■  <88> 

zo 


Hence  to  obtain  V  ,  the  diffuse  scatter  p ,  (0_)  is  scaled  by  the  factor 
s  d  Z 

A'/(2 rr  v3  ),  and  the  abscissa  is  converted  from  0_  to  k  sin  0„.  Because 
zo  z  Z 

of  isotropy  of  statistics,  the  relation  is  independent  of  0^  and  requires  a 
measurement  of  scattered  power  for  Os  0  <  tt/2  for  any  0  .  Even  though 

b  J 

V  (v)  has  non-zero  values  for  v>k  (where  v  =  k  is  the  limit  of  measure- 
s 

ment),  we  must  make  some  assumption  on  the  character  of  V  (v>k).  Either 

s 

V  (v>k)  is  set  to  zero,  or  it  is  set  by  some  trend  (e.g.  Gaussian)  which  is 
s 

established  for  v<  k. 

In  appendix  III,  the  one-dimensional  power  spectral  density  P 

s 

is  related  to  V  by 
s 


00 


(89) 


PjO  =  J 

K 


-  -2  - -  V  (v)  dv 

Jv3  -  K3  S 


where  tc  =  2  tt  v  .  (The  only  restriction  on  Vg  is  that  it  be  Bessel-transform¬ 
able.)  Unfortunately  the  transform  in(89)  does  not  possess  an  inverse  which 
is  a  single- variable  integral,  and  tables  for  this  transform  have  not  been 

generated.  Hence  the  most  straightforward  way  to  go  from  V  to  P  is  to 

s  s 

Fourier- Bessel  transform  V  ,  and  then  Fourier- cosine  transform  the  result. 

s 

To  go  from  P  to  V  ,  the  procedure  is  reversed.  There  is  one  caution  in 
s  s 

getting  V  from  P  :  that  the  Fourier- cosine  transform  of  P  must  be  Bessel- 
s  s  s 

transformable.  This  immediately  eliminates  the  Elrac-delta  function  (for  P  ) 

s 

which  possesses  a  non- convergent  Fourier- cosine  transform.  However,  if 

V  is  a  delta  function  of  area  2c: 
s 


Vg(v)  --  c  r  «(v- Vj )  +  6  ( V  +  Vj )  ] 


then 


P  (K)  = 
s 


2  c  v. 


•J  V3  - 


1*1  s  Vj 

I  *  I  >  v! 


(90) 


(91) 


This  says  that  information  stored  as  a  delta  function  at  v  =  +  v,  in  V  affects 

-  1  s 

the  character  of  P  (  k  )  for  1  k  I  s  v, ,  but  not  for  I  K  I  >  v, . 

s  1  '  1  ’  1  l 

Relations  (90)  and  (91)  show  how  to  choose  a  suitable  Nyquist 

frequency  when  P  (v)  is  measured  by  a  surface  tracing  method  (e.g.  on  a 
s 

proficorder) .  The  maximum  value  on  the  abscissa  of  V  is  v  =  k  (since 

s 

|  09  |  <  tt/2).  In  the  abscissa  of  P  ,  that  maximum  corresponds  to  u=  2nk  = 
-I  S 

X  .  Hence  for  a  scatter  measurement  for  all  0^  (rather  than  Just  0^0), 
a  suitable  Nyquist  frequency  is 

vny=2X_1*  (92) 

(Namely  the  sampling  period  would  be  X/4  or  shorter.  )  Alternatively,  a 

determination  of  P  via  a  scatter  measurement  can  be  made  for  frequency  v 
s 

greater  than  physically  possible  with  the  most  sensitive  surface  tracing 
devices.  For  example,  generally  on  the  best  proficorder  surface  trace,  a 
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minimum  reasonable  sampling  period  is  about  0.1  mil  (which  corresponds 
to  a  Nyquist  frequency  of  0,.  196  p,  *  ).  However,  for  a  Nd-YAG  laser  beam 
(X=  1.06  |j,)  scattered  from  the  same  surface  (at  normal  incidence),  the 
Nyquist  frequency  is  0.943  y,  *.  (Of  course,  inequality  (83)  must  still 
be  satisfied.) 

To  see  how  to  use  equations  (88)  and  (89),  we  look  at  a  concrete 
example  by  returning  to  the  Gaussian  model  fdr  the  one-dimensional  power 


spectrum: 


P  (  u  )  =  a2  t  Jr T  bxp  [  -(  vtt  t  )2  ]  . 
s  so  o 


(A  one-surface  model,  n  =  1 ,  is  sufficient  for  explanatory  purposes  and  can 

s 

easily  be  extended  to  higher  orders,  )  Hence  we  find 

rr  r  2  a  3  v  3 

pd(v)  =  - d---— 2-°  exp.  [-(tqv/2)3  ]  .  -  ,  (94) 

‘  I 

If  the  Gaussian  model  is  appropriate  for  P  ,  the  surface  parameters  a  and 

s  s 

t  can  be  found  from  a  quick  examination  of  the  scatter  pattern.  Initially, 
o  ' 

the  peak  of  the  diffuse  component  p,(0)  gives 

,,  d 

O  P H(0)  A'  1  ! 

(  cr  t  )  =  -  •  05) 

SO  4  k*3  TT  ‘  I 

(The  A'  is  found  in  the  analysis  of  the  coherent  component.)  Secondly, 

suppose  the  e-folding  angle  of  p,(0o)  is  0  o  ,  and  that  09'  <  tt  /4.1 
/  a  2  *0  40 

Then  to  a  close  approximation 


To  k  0, 


and  of  course  (96)  in  (95)  gives  a  .  We  can  extend  this  procedure  to 

!  i 

higher-order  composite  surfaces  by  examining  each  component  of  scatter 
individually  (provided  the  components  are  distinguishable).  The  e-folding 

i 

angles  of  the  coherent  and  diffuse  components  (assuming(93)  is  the  power 
spectrum)  are  depicted  in  the  hypothetical  scatter  plot  of  figure  14. 


I 


Figure  14.  Hypothetical  scatter  at  normal  incidence  from  a  slightly  rough 
surface  with  Gaussian  p.s.d.f.  and  p.d.f.  showing  relation 
between  e-folding  widths  and  peaks  of  components . 


Figure  15.  Hypothetical  backscatter  (otherwise  like  figure  14). 
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2 .  Direct  Backscatter. 

There  is,  however,  another  measurement  which  can  be  made  to 

give  p  ,  t  ,  and  a  and  in  some  cases  the  surface  power  spectral  density, 
o  o  s 

If  instead  of  the  scatter  for  vertical  incidence  being  measured,  suppose  that 
the  direct  backscatter  from  the  surface (0=02=-9^/  83=  0  )  is  a  more  conven¬ 
ient  measurement.  The  expression  analogous  to  (84)  is 


2  2  « 

/  n\  T  A'v  l  L  4  k  secfl  i2 vx  . 

PB(  6)  =  exp  | —  J  +  - -  j  J  Bs(t)  e  dx  dy . 


Unfortunately  (97)  does  not  provide  us  a  means  of  finding  a  general  B  (t) 

s 

from  p  (8) .  If  however,  the  covariance  can  be  decoupled  as 

D 

B  (t)  =  a2  C  (x)  C  (y)  ,  (98) 

s  s  s !  si 

then  the  covariance  can  be  found  from  a  measurement  of  diffuse  power,  Prf(9) 
(the  second  term  in  (97)  ).  If  we  define  the  Fourier-cosine  transform 


I(x)  =  f  (cos  9)  Pj(9)  cos  (  2  v  x  )  dv 


then  C  (x)  =  I  (x)  /  I  (0) 

S1 


r  (T)  =  A'  I  (x)  I(y) 

s  T  2 

4  tt  k  J  I  (x)  dx 


(100) 


(101) 


From  the  covariance,  either  the  one  or  two-dimensional  power  spectral  den¬ 
sities  can  be  found  (as  Fourier-cosine  or  Fourier- Bessel  transforms  of  B  (r), 

s 

respectively).  As  in  the  case  of  scatter  for  normal  incidence,  the  constant 

A'  is  found  from  the  coherent  component  of  scatter,  p  (0)  (first  term  in  (97)  ). 

c 

If  the  e-folding  angle  of  the  coherent  component,  6  ,  is  much  less  than  rr/2, 

o 

then 


A'  = 


(102) 


Moreover,  the  e-folding  length  of  the  OTF  is 
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2 

and  C  is  found  from  (87).  If  in  addition  the  Gaussian  model  (93)  for  P 
n  s 

is  assumed,  and  the  e-folding  angle  of  the  diffuse  component  of  scatter, 

9'  ,  is  less  than  about  tt/4,  then  to  a  close  approximation: 

t  =  (  k  9  '  f 1  (104) 

o  o 


and 


(cr  t  )‘ 
s  o 


P.(0)  A' 

d 

4  k5*  TT 


(105) 


The  e-folding  angles  and  peaks  of  the  diffuse  and  coherent  components  are 

depicted  in  figure  15.  From  a  direct  comparison  of  figures  14  and  15,  we 

see  that  backscatter  p  (0)  is  exactly  the  same  as  the  direct  scatter  with  its 

B 

abscissa  compressed  by  a  factor  of  two,  p(  2  9^  ).  This  conclusion  does 
not  apply  in  general  but  is  at  least  valid  for  a  Gaussian  power  spectral 
density.  Extension  to  higher  order  surface  composites  is  straightforward 
provided  the  individual  components  are  distinguishable. 
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Chapter  IV 

EXPERIMENTAL  RESULTS 

This  chapter  presents  measurements  of  scattering  of  a  CC>2  laser 
beam  from  the  aluminum  and  titanium  surfaces  which  were  analyzed  on  the  pro- 
ficorder.  The  turbulence  is  assumed  negligible,  and  the  surfaces  have  a  zero 
mean;  hence  the  comparison  of  experimental  results  and  theoretical  computer¬ 
generated  curves  will  serve  as  a  verification  of  the  rough  surface  and  profile 
aspects  of  the  theory.  In  addition,  this  checks  the  accuracy  of  the  composite- 
surface  model.  There  are  four  patterns  of  scatter  in  the  plane  of  incidence 
(  9=0)  presented:  scatter  from  titanium  for  9  =27°  and  3.5°,  and  from  alum- 

X 

inum  for  18°  and  3.5°  .  (The  beam  was  incident  upon  the  location  where  pro- 

ficorder  traces  were  made.)  The  data  for  the  near-normal  incidence  (0  ^3.5°) 

are  also  used  to  derive  the  surface  Pg  via  the  procedure  of  the  "inverse  problem." 

(Achieving  0^0  is  physically  impossible  in  the  present  experimental  set-up.) 

As  a  further  check  of  the  rough  surface  model,  the  power  spectrum  PS  (super-s 

®  t 

to  indicate  originated  from  a  scattering  measurement)  is  compared  with  P 

s 

(super-t  for  proficorder  trace) . 

A.  EXPERIMENTAL  ARRANGEMENT 

This  section  outlines  the  apparatus  used  in  making  the  scatter 
measurements.  (See  figure  16.)  The  CC>2  laser  beam  passes  through  a  beam 
splitter  (for  monitoring  laser  power  on  the  thermopile) ,  through  a  chopper 
wheel  (for  lock-in  amplification) ,  and  is  reflected  to  the  scattering  surface 
mounted  on  the  target  table.  All  components  are  adjusted  until  the  laser  beam, 
the  normal  direction  to  the  surface,  and  the  detector  are  coplanar.  The  max¬ 
imum  distance  from  target  to  detector  is  50  cm;  and  a  lens  at  the  detector 
increases  signal  strength.  The  detector  is  mounted  on  an  optical  bench  which 
is  rotated  manually  about  the  scattering  surface.  Below  the  bench  is  a  pointer 
which  indicates  the  scattering  angle  9  with  accuracy  better  than  +  £  degree. 

The  major  system  components — the  laser  and  its  monitor,  the  detector  and  its 
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Figure  16a.  Experimental  set-up  for  measurement  of  scattered  power  versus 
$2  when  63  =  0  and  G1  is  set. 


lens,  and  the  electronics  —  are  now  described  individually. 

The  CC^  laser  built  by  the  author  has  a  CW  output  power  of 
40  watts  (11%  efficiency) .  Beam  diameter  at  the  output  mirror  is  approximately 
1.2  cm  (as  measured  at  e-folding  points  of  the  power) ,  and  divergence  is 
0.4  mrad.  Three  gases  —  CC>2  at  a  pressure  of  1.4  torr,  Ng  at  1.6  torr,  and 
He  at  7  torr  —  are  flowed  through  the  one  inch  inside  diameter  pyrex  tube. 

The  gold-coated  pyrex  mirror  (diameter  3.81  cm,  radius  of  curvature  10m) 
and  the  germanium  output  mirror  (diameter  2.54  cm,  flat,  80%  reflection) 
are  at  a  separation  of  130  cm.  Irtian  II  windows  at  the  Brewster  angle  give 
a  vertical  polarization  at  the  target  surface.  The  active  discharge  region  is 
86  cm  in  length  and  is  cooled  by  a  concentric  water  jacket.  The  output  of  the 
power  supply  (operated  at  60ma,  5500v)  is  filtered  to  give  0.05%  voltage 
ripple  as  measured  at  the  discharge  tube.  The  major  source  of  short-term 
laser  noise  is  the  uneven  pumping  of  the  Welch  fore-pump.  Long-term  varia¬ 
tions  are  a  result  of  the  heating  of  the  mirrors  which  caused  the  laser  to 
tune  and  de-tune  itself.  Generally  the  long-term  fluctuations  stabilize  after 
two  hours  of  operation.  Since  the  precise  laser  power  is  not  certain,  the  power 
is  constantly  monitored  by  reflecting  part  of  the  beam  from  an  Irtran  II  flat  to 
a  thermopile  (Eppley  8-junction  Bismuth-Silver) .  All  measurements  of  scattered 
power  are  normalized  to  the  laser  power  as  measured  at  the  thermopile. 

The  detector  is  a  Santa  Barbara  Research  Corporation  gold-doped 
germanium  photoconductive  device.  All  measurements  of  power  were  of 
sufficiently  low  level  that  the  detector  always  operated  in  the  linear  regime 
of  its  bias  circuitry.  A  thin  uncoated  germanium  wafer  (transmitting  10%) 
protects  the  detector  while  the  coherent  component  is  measured.  Also  a 
germanium  meniscus  lens  of  diameter  2. 5cm  and  focal  length  2.54  cm  focuses 
the  radiation  onto  the  detector.  Over  the  face  of  the  lens  is  a  diaphragm  which 
can  vary  the  f/stop  from  f/1  to  f/10  continuously  for  a  more  selective  sampling 
of  laser  power.  The  lens  serves  four  functions:  (a)  it  increases  the  signal 
level  at  the  detector;  (b)  it  increases  the  signal  to  noise  ratio  by  filtering  the 
extraneous  near-infrared  and  visible  light;  (c)  it  has  the  effect  of  reducing  the 
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distance  *o  the  far-field;  and  (d)  it  performs  the  surface  ensemble  average  as 
required  by  the  theory. 

The  system  electronics  consist  of  two  oscilloscopes,  a  lock-in 
amplifier,  two  pre-amplifiers,  and  a  chopper  wheel.  A  Tektronix  503  oscil¬ 
loscope  monitors  the  thermopile  voltage.  (A  micro-voltmeter  could  be  used  as 
well.)  The  signal  from  the  detector  is  fed  to  a  PAR  211  pre-amplifier  and  then 
to  a  PAR  120  lock-in  amplifier.  The  pre-amplifier  matches  the  detector  to  the 
lock-in  amplifier  without  loading  the  detector  bias  circuit,  and  it  allows  vari¬ 
able  amplification  from  zero  to  60  db.  The  output  of  the  pre-amplifier  is 
monitored  on  an  oscilloscope  for  visual  inspection  of  the  signal.  For  a 
reference  signal,  a  zener  diode  (detector) views  a  pilot  lamp  through  the 
chopper  wheel  (frequency  330  Hz).  That  signal  is  amplified  in  a  Tektronix 
1A7A  plug-in  and  fed  as  a  "selective  external"  input  to  the  lock-in  ampli¬ 
fier.  The  lock-in  amplifier  is  needed  to  measure  the  diffuse  component  of 
scatter  since  only  the  coherent  component  is  easily  measured  on  the  oscil¬ 
loscope.  The  lock-in  amplifier  has  a  time  constant  set  to  one  second  to 
average  the  laser  short-term  fluctuations.  This  value  corresponds  to  the 
one  second  time  constant  of  the  thermopile  monitoring  the  beam  power. 

B.  MEASUREMENTS 
1.  Preliminary  Comments. 

Measurements  were  made  in  the  near- field  of  the  surfaces  because 
of  the  physical  difficulties  (e.g.  alignment)  of  positioning  the  detector  in 
the  far-field  and  because  of  the  enhanced  signal  level  in  the  near- field. 

(The  far-field  is  for  ranges  Rq  »  50m  for  a  beam  diameter  at  the  surface  of 
1  cm.)  Good  agreement  for  this  set-up  is  expected  for  the  diffuse  component 
of  scatter,  but  not  for  the  coherent  component.  The  coherent  scatter  is  not 
yet  fully  developed  when  the  target-to-detector  distance  is  still  significant 
with  respect  to  beam  diameter.  Most  important,  in  all  measurements  of  scatter 
(figures  17-20)  ,  the  coherent  component  measured  is  equivalent  to  the  scatter 
from  a  perfectly  smooth  surface  replacing  the  target  surface.  (On  each  figure 
a  Gaussian  is  fitted  through  coherent  data  to  show  where  the  coherent  component 
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stops  and  the  diffuse  starts.)  The  measurements  are  not  needed  to  verify 
the  coherent  scatter  since  in  part  4  ,  section  F  of  chapterlll  we  see  the 
reasonability  of  the  theoretical  form  of  coherent  scatter. 

The  measurement  is  of  course  power  versus  02  for  0^  and  0^  con¬ 
stants.  However,  the  theory  gives  power  flux  (density).  A  comparison  of 
data  with  theory  is  still  legitimate  provided  the  size  of  the  detector  aperture 
does  not  exceed  large  scale  spatial  power  variations.  Small  scale  variations 
will  be  integrated  since  that  average  corresponds  to  the  large  ensemble  aver¬ 
age  taken  in  the  theoretical  derivation.  Hence  power  density  curves  of  theory 
are  compared  with  experimental  power  curves  by  adjusting  the  ordinate  of  one 
with  respect  to  the  other.  No  such  adjustment  of  abscissa  is  necessary 
since  peaks  of  coherent  components  indicate  a  reference  angle. 

2 .  Titanium  Scatter. 

The  data  and  theoretical  curves  for  scattering  from  titanium  where 
01=27°  are  in  figure  17  and  where  0^3.5°  in  figure  18.  General  overall 
agreement  is  obvious,  especially  for  |  62“ei  I  ^  n/4  ;  however  three  dis¬ 
crepancies  will  be  accounted  for  in  the  order  of  increasing  ©2  from  specular. 
First,  a  slight  ripple  is  noted  at  @2»  5°  and  57°  in  figure  17  and  at  ©2  »  37° 
in  figure  18.  We  can  see  the  source  of  this  seeming  anomaly  in  the  power 

spectrum  P*  (figures  6  and  8)  at  v  «0.05n  *.  Irregularities  of  this  sort 
s 

can  not  be  accounted  for  by  the  Gaussian  model,  but  they  serve  the  more 
important  role  of  providing  a  target  signature.  A  signature  is  a  characteristic 
of  scatter  from  a  surface  which  is  not  seen  for  other  surfaces. 

Second,  for  values  |  0  —0  |  >  40°  in  the  roughness  component  of 

X  M 

scatter,  the  data  are  higher  than  theory  predicts .  The  discrepancy  is  not 
attributable  to  any  theoretical  error  because  no  assumptions  are  violated 
until  |  01 1  or  I e  2  I  aPProaches  90°  .  Moreover  such  an  error  does  not  appear 
in  the  aluminum  measurements  (figure  20).  Rather  the  explanation  is  that 
the  stylus  of  the  proficorder  can  not  make  a  perfect  trace  of  the  surface  but 
will  slightly  indent  the  material.  This  problem  will  become  worse  with 
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Figure  17.  Measured  scattered  power  and  theoretical  power  flux  for 

titanium,  8^  =  27°,  8^  =  0.  (The  vertical  scale,  in  db  normal¬ 
ized  to  backscattered  power  at  vertical  incidence,  in  figures 
17-20  is  fcr  the  theory;  ordinate  for  data  is  arbitrary.) 
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increasing  wavenumber  in  the,  power  spectrum.  Hence  the  measured  values 
of  Ps(v)  for  v  ^vjsjy//^  wi^  ^e  lower  than  is  actually  the  case.  This  problem 
is  not  so  bad  for  aluminum  s  ince  the  aluminum  potoer  spectrum  is  at  least  a 

.  i  ,  < 

factor  of  two  lower  in  this  large  wavenumber  region. (See  figures  8  and  9.) 

Finally  ,  for  scatter  near  the  surface  edge  ,  the,  data  fall  off  more 

/ 

rapidly  than  is  predicted.  Recall  that  the  theory  does  not  apply  for  grazing 
,  incidence  ( 1 0  |  «  tt/2)  or  grazing  scatter  (  |@  |  «  rr/2  ).  Appendix  I  analyzes 

i  ^  ;  »  ^ 

the  grazing  problem.  •' 


3 .  Aluminum  Scatter. 

The  data  and  theoretical  curves  for  scatter  from  the  (corroded)  1 

i 

aluminum  where  0^  =  18°  are  given  in  figure  19  and  where  0^  =  3.5°  in  figure 
20.  The  scatter  of  figure  19  is  especially  provided  for  the  waviness  component 

f  f  .  ' 

,of  scatter;  whereas  figure  20  depicts  the  excellent  agreement  in  the  roughness 
component.  Contrary  to  titanium  measurements,  the  waviness  component  does 

'  I  '  I  ; 

not  exhibit  good  agreement.  Initially  the  modeled  value's  for  a  and  t  are 

'  '  '  '  si  °1 

suspect.  From  the  data  of  figure  20  and  the  procedure  of  the  inverse,  problem 

(assuming  a  Gaussian  P„  ) ,  we  find  o  =  0.044  u  (rather  than  the  modeled 

■  S1  S1  ■ 
0.056  m,)  and  tq  =  37.0  p,  (rather  than  the  24. 9  n,  used) .  (The  vailueS  ^re  found 

without  a  knowledge  of  the  constant  A'  by  assuming  the  values  of  a  and.T 

.  s2  °2 

are  correct  and  then  noting  the  ratio  of  the  peaks  of  the  individual  diffuse 

components  pd^(0)  /  p^  (0)  equals  (as^T0  )  / . (crs  tq  )  .)  If  the  corrected 

values  of  a-  and  x_  were  used,  the  width  of  the  waviness  component  would 
S1  °1  r  ,  ■ 

be  smaller  and  its  peak  higher. 

,  f 

In  the  sectiop  presenting  surface  measurements  (chapter  Ilf,  section 

>  i 

A;  and  figure  7)  we  saw  how  traces  at  different  places  on  the  surface  give 
quite  different  power  spectra  .  The  non-stationarity  does  not  significantly 
affect  the  scattering  behavior  provided  the  beam  covers  the  largest  scale 
surface  variations  or  provided  the  beam  size  is  so  small  that  the  surface 

i  '  i 

under  the  beam  is  actually'  stationary.  When  scattering  from  an  airplane 
surface  at  a  large  range,  the  beam  diameter  will  be  large  enough  (about  one 

I  *  . 

meter  or  greater)  to  average  the  non-stationarity  so  that  the  scattering  will 
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Figure  19.  Measured  scattered  power  and  theoretical  power  flux  for  aluminum, 
ei  =  18°'  03  =  O* 
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appear  stationary.  In  the  laboratory,  the  1  cm  beam  diameter  exposes  the 
non-stationarity  of  the  surface. 


4.  Titanium  Inverse  Measurements. 


From  the  data  of  figure  18,  the  two-dimensional  power  spectrum 

s 

V  and  the  one-dimensional  power  spectrum  P  are  determined.  Since  only 
s  s 

the  diffuse  scatter  is  needed,  the  first  three  data  points  in  p^  are  estimated. 
As  demonstrated  by  equation  (88) ,  Vg  is  found  by  a  scaling  of  p^ .  (That  spec¬ 
trum  is  plotted  in  figure  21.)  After  consecutive  Bessel  and  Fourier  transforms, 
PS  is  found.  In  figure  22  are  PS  and  P1  .  Since  the  proper  scaling  of  p,  is 

S  S  S  Q 

not  known  (due  to  the  inability  to  determine  p  in  the  near- field) ,  there  is  an 

s  c  s 

arbitrary  constant  which  multiplies  Pg .  Hence  the  ordinate  of  Pg  is  adjusted 

until  a  good  match  is  obtained  with  P*  .  The  scale  factor  would  be  known  if 

2  S 
a  were  measured  by  some  other  method.  In  the  laboratory  this  poses  no 
s 

problem  since  inexpensive  mean- height  indicators  are  generally  available. 

If  instead,  scatter  from  an  airborne  target  is  measured,  the  receiver  would 
be  in  the  far- field,  and  once  again  the  variance  would  be  determined. 


TITANIUM 


Figure  21.  Titanium  two-dimensional  power  spectral  density  V  versus 

S  4 

wavenumber  k  .  (Ordinate  has  units  proportional  to  u  .) 
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t  s 

Figure  22.  Titanium  P  (surface  height  p.s.d.f.  from  tracing)  and  P 

s  s 

(surface  p.s.d.f.  from  scattering).  P  is  found  from  the  data 

s 

of  figure  18  and  the  procedure  of  the  inverse  problem. 
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From  a  comparison  of  P  and  P  ,  we  note  the  error  in  the  form  of 

£  S  S 

P  (v)  is  for  frequency  v  approaching  Nyquist.  As  noted  earlier,  the  discre- 
pancy  is  due  to  the  inability  of  a  proficorder  to  trace  small-scale  surface  vari¬ 
ations  of  significant  height  accurately.  But  when  the  power  spectrum  is  mea¬ 
sured  by  scattering,  the  surface  is  not  deformed;  hence  the  number  of  samples 
(all  points  under  the  beam)can  be  large,  and  stationarity  can  be  easily  verified. 

5.  Aluminum  Inverse  Measurements. 

The  aluminum  two  and  one-dimensional  power  spectra  (figures  23 
and  24  respectively)  are  obtained  from  the  data  of  figure  20  by  the  same  steps 

g 

as  for  titanium.  The  plot  of  P  better  demonstrates  the  error  in  the  waviness 

s  t 

structure.  If  the  second  through  fourth  points  in  Pg  had  been  used  to  estab¬ 
lish  the  waviness  model,  the  agreement  between  scattering  theory  and  data 
vwuld  have  been  excellent.  From  an  inspection  of  the  individual  traces  of 


*  (  /i"'  ) 

Figure  23.  Aluminum  two  dimensional  power  spectral  density  V  versus 

s 

wavenumber  K  .  (Ordinate  has  units  proportional  to  y,4.) 


Figure  24. 


P*  in  figure  7,  we  see  that  a  single  plot  (number  3;  13%  of  all  the  data) 
s 

caused  the  error  in  waviness  component.  This  indicates  that  because  of  the 
non-stationarity  of  aluminum,  a  very  large  number  of  samples  must  be  taken 
to  get  an  average  power  spectrum. 

6.  Conclusions. 

Taking  a  two-surface  Gaussian  model  for  the  surface  power  spectrum 

g 

is  justified  by  the  form  of  P  for  titanium  (figure  22)  and  for  aluminum  (figure 

s 

24).  In  those  plots,  the  waviness  and  roughness  components  are  even  more 

distinct  than  in  the  plots  of  P*.  The  improvement  is  of  course  due  to  the  large 

number  of  samples  the  beam  takes  when  it  scatters  from  an  area  on  the  surface. 

To  get  a  smoother  plot  of  P*  ,  the  number  of  data  points  should  be  increased 

s 

to  at  least  8000  for  aluminum  and  2000  for  titanium  (because  of  the  more- 
nearly  stationary  character  of  the  latter) .  But  still  the  problem  of  deformation 
of  the  surface  by  the  stylus  cannot  be  avoided,  and  the  increased  sampling 
will  not  help  the  titanium  agreement.  To  check  for  surface  stationarity ,  the 

3 

beam  size  might  be  varied  and  Pg  found  for  each  size.  Conceivably  there  is 
a  size  beyond  which  the  diffuse  pattern  does  not  change  further.  For  that  size 
the  surface  appears  stationary,  and  an  average  power  spectrum  is  the  proper 
one  for  use  in  modeling . 

From  the  scattering  measurements  already  done,  we  can  visualize 
an  optimal  system  for  the  determination  of  rough  surface  (isotropic,  homogen¬ 
eous)  statistics.  First  the  probability  density  is  found  on  a  multi-channel 
analyzer  (MCA)  by  evaluation  of  the  proficorder  output.  If  the  density  is  nearly 

Gaussian  and  isotropic,  then  the  theoretical  derivation  of  the  inverse  problem 

2 

is  appropriate.  Since  a  proficorder  is  being  used,  the  surface  variance,  a  , 

S 

is  easily  measured  —  either  from  the  meter  on  the  proficorder  or  as  the  second 
central  moment  of  the  probability  density.  Next  the  appropriate  laser  is  chosen 
for  scattering,  depending  on  what  Nyquist  frequency  (\  *)  is  desired.  The 
beam  is  to  have  an  incident  angle  as  near  normal  as  possible.  The  scattered 
pattern  from  specular  to  the  edge  is  measured  and  fed  to  a  computer  program 
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s  s 

which  scales  and  transforms  to  get  P  .  Final  scaling  of  P  is  done  with  the 
2  s  s 

a  found  by  another  method.  Finally,  from  the  knowledge  of  P  ,  one  might 
s  s 

filter  the  proficorder  output  into  waviness  and  roughness  components  for 

individual  analysis  on  the  MCA.  The  advantage  of  this  system  over  recording 

the  trace  and  analog-to-digital  converting  is  in  time  saving.  An  excessive 

number  of  data  points  must  be  used  to  properly  establish  the  power  spectra 

via  tracing,  but  only  about  90  to  180  data  samples  are  taken  in  scattering. 

If  the  cost  of  the  MCA  and  proficorder  is  prohibitive,  then  the  investigator 

might  establish  the  probability  density  by  theoretical  considerations  (central 

2 

limit  theorem)  and  measure  a  on  a  simple  (inexpensive)  roughness-measuring 

s 

device.  In  any  case,  the  computer  is  only  used  for  data  reduction,  and  the 
problem  of  analog-to-digital  conversion  is  avoided. 


Chapter  V 
CONCLUSIONS 

With  the  increasing  development  and  use  of  new  airplane  surfaces 
(such  as  Graphite/Epoxy  and  Boron/Epoxy) ,  measurements  of  scattering  and 
the  establishment  of  new  models  will  need  to  be  effected  for  the  n*w  surfaces. 
Also  one  Surface  presently  in  use,  the  painted  one,  needs  to  be  modeled  by 
persons  who  have  the  material  available  to  them.  Since  additional  information 
is  available  on  the  surface  when  the  wavelength  is  changed,  scattering 
measurements  should  be  made  at  wavelengths  shorter  than  the  10. 6  y,  used  in 
this  analysis.  (As  pointed  out  in  the  text,  the  aluminum  and  titanium  surfaces 
are  moderately  rough  at  1.0  M-;  hence  the  theory  should  still  give  good  agree¬ 
ment  there.)  For  all  surfaces  examined,  the  power  spectra  and  scattering  curves 
should  be  carefully  inspected  for  unique  characteristics  which  will  be  useful 
in  target  identification.  Such  a  characteristic  was  found  for  titanium  as  a 
"ripple"  in  the  power  spectrum. 

The  theory  (in  the  form  of  equation  (66)  )  is  now  ready  for  the  gen¬ 
eration  on  a  digital  computer  of  simulated  target  images.  The  target  would 
be  described  by  grid-points  in  a  three-dimensional  space,  where  the  spacing 
of  grid  lines  is  less  than  the  beam  diameter.  A  quadric  surface  v\ould 
approximate  the  deterministic  surface  at  each  grid  crossing,  and  the  scatter 
to  the  detector  at  some  bistatic  or  monostatic  angle  would  be  calculated.  The 
assemblage  of  the  returns  constitutes  the  image.  Then  the  techniques  of 
pattern  recognition  would  be  applied  to  help  establish  a  target-image  relation¬ 
ship.  At  the  same  time  laboratory  measurements  of  target  images  will  further 
establish  the  relation. 

Turbulence  was  carefully  injected  into  the  scattering  problem  because 
a  disturbance  of  the  incident  plane  wave  character  will  obviously  affect  the 
scattering.  However  ,  scintillation  of  the  beam  on  the  return  path  is  not  an 
integral  part  of  the  scattering  and  hence  must  be  treated  separately.  Once 
again  the  optical  transfer  function  would  be  used  but  would  be  applied  to 
the  total  target  image  in  the  image  coordinate  system.  The  optical  transfer 
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function  derived  in  the  text  is  not  limited  to  use  in  this  scattering  problem. 

The  atmospheric  OTF  is  important  in  all  problems  of  propagating  a  coherent 
beam  through  a  turbulent  medium.  The  result  will  find  uses  in  communications 
and  radar. 

Finally,  the  reader  should  be  especially  aware  of  the  transfer  function 
formulation  which  allows  the  investigator  to  vary  some  aspects  of  the  problem 
without  disturbing  the  rest.  The  solution  was  kept  as  general  as  possible 
before  applying  it  to  laser  scattering  from  rough  metallic  surfaces  so  that 
the  solution  can  be  easily  applied  to  other  scattering  problems,  including 
underwater  scattering  from  targets  and  scattering  from  the  ocean  surface. 
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APPENDIX  I 

Verification  of  Inequalities  (11)  and  (18) 


This  appendix  presents  more  explicit  inequalities  to  replace  (11) 
and  (18)  based  on  the  chosen  field  profile  and  OTF  of  chapter  III.  Neither 
inequality  can  be  strictly  satisfied  since  U  is  a  stochastic  quantity  which  can 
go  to  zero.  However,  the  replacement  yields  inequalities  not  so  precise  but 
easier  to  use. 

For  an  isotropic  homogeneous  random  variable  with  Gaussian 
autocorrelation,  the  standard  deviation  (o')  of  the  derivative  of  the  variable 
is  given  by 

o'  =  o  /  T  (106) 

o 

where  a  is  the  standard  deviation  of  the  variable,  and  T  is  the  e-folding 

o 

distance  of  the  autocorrelation.  Also  the  following  replacement  is  used: 

-y*  =  (w+icp)r+~*  (107) 

where  w  =  In  |z|  and  cp  =  arg  (Z)  (108) 


represent  the  random  nature  of  the  medium,  and  the  sub-r  indicates  the  x,y 
or  z  partial  derivative.  Using  (107) ,  each  inequality  is  split  into  a  "deter¬ 
ministic"  and  a  "statistical"  inequality.  From  (106),  the  approximations 


|(w  +  i  <p)x  |  *  (cos  9^  vT  /  PQ  , 

|  (w  +  i  cp)y  I  <  J2  /  pQ 

|(w+i«p)2|  ^  I  sin  |  J  2  /  pQ  (109) 

— 1/8 

are  made.  (Note  the  variance  of  (w  +  i  9)  is  Bu(0) ,  and  To=  0.945  Lq  [B^ (0)  ]  . ) 

In  a  similar  fashion: 


^  (  cos  6  j)  2/X  , 
*  2/1  , 
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and 


e 


^  •  sin  0  |  2/1. 


(110) 


Using  (107),  (109),  and  (110),  the  first  inequality  (11)  decouples 
into  a  "deterministic"  inequality 


2 

kJt 


tan  6  |  «  1 


(111) 


and  a  "statistical"  inequality 


kf ■'l“n»1  I  <<!• 


(112) 


Inequality  (18)  weakens  when  9^  =  0  and  (specular  scatter 

near  grazing).  For  these  angles,  (18)  becomes  ,  1 


h  se°2  ei  ^ 1 


and 


1  sec2  e  «  1 . 

k  Pc  1  ' 


(113) 

(114) 


In  summary,  the  four  inequalities  (111-114)  depend  only  on  the  Wave- 

:  .  I 

length,  the  incident  angle,  and  the  e-folding  lengths  of  the  OTF  and  the  inci¬ 
dent  field.  The  inequalities  say  the  work  of  this  paper  is  not  applicable  to 

f 

the  study  of  scattering  for  grazing  incidence.  Finally,  use  of  the  weak  bounds 
of  (109)  and  (110)  is  justified  in  that  the  inequalities  are  indeed  ejasiiy  satis¬ 
fied.  If  the  inequalities  had  turned  out  to  be  marginal,  then  a  more  precise 

j  '  i 

replacement  would  have  been  necessary.  ,  , 
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Integration  for  the  Wave  Correlation  Function 

,  '  i 

This  section  briefly  outlines  the  steps  taken  in  arriving  at  the 
approximate  form  (55)  for  the  wave  correlation  function.  Initially  the 
integral  (52)  which  defines  the  correlation  is  rewritten  as 


h_L~ 

mo 


C^(p')~  J  JoU'p'.)  (1+  h'2)  "11/6  h'  dn' 


(115) 


(where  p'  =  p/LQ  ,  and  =  kLq  )  to  reduce  the  number  of  parameters  from 

two  (  k  and  L  )  to  one  (k  L  ).  Equation(l  15)  was  numerically  integrated 
m  o  mo  ■ 

using  the  rectangular  rule  on  a  digital  computer.  It  was  normalized  such  that 

C  (p')  at  some  small  p'  (less  than  1/(k  L  )  but  not  zero)  matches  the 
u  K  mo 

value  for  (  p')  as  given  by  (54) .  Next  the  function 


y  =  -  _£n  Cu  (p ',) 


(116) 


was  plotted  on  log-log  paper  in  order  to  determine  the  parametric  dependence 
and  possible  functional  dependence  (perhaps  C  ~  exp  (-b  p'  a) ,  b  and  a 

i  ti 

constants) .  That  plot  (of  which  a  significant  portion  is  in  figure  25)  demon¬ 
strates  the  exponent/a,  decreases  from  2.0  to  1.5  as  p'  increases  from 


-4  -1 

10  to  10  , 


However,  a  Gaussian  form  for  C  is  desired;  hence  a  suitable 
’  u 


method  of  approximating  a* 2  is  presented.  Since  form  (54)  is  essentially 

valid  for  p'  <  1/  (h  L  ),  the  following  curve  (117)  is  also  plotted  (the  dashed 

mo 

lihe)  as  an  upper  boundary  for  the  validity  of  (54): 


/  /\  _  ^r/  * \ ”  1/3  ^  i  (  »v2 

yi'p  )  -4  l(p  )  -  j  J  (p  )  . 


Then 


yi(p')  «  0.695  (p')8/5 


10“5<p'<  10_1 


(117) 


(118) 


is  found  as  a  graphical  approximation  to  (117).  Fortunately  for  a  constant 
Bu(0)  and  the  OTF  Cq  varying  from  0, 9 to  0.1 ,  the  function  y  =  -£n  (C^) 
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Figure  25.  Plot  used  in  determining  the  approximate  Gaussian  form  for 
the  OTF.  (Symbols  correspond  to  computer  generated  data. 
See  appendix  II.) 
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varies  only  slightly  more  than  one  decade.  (For  a  given  Cq,  say  0.9  or 
0.1,  the  ordinate  is 


yo  ~ 


I  ^  Cp  | 

Bu») 


(U9) 


Then  note  An  (0.1)/  An  (0.9)  =  22.)  Over  this  decade  (termed  the  "operating 
region") ,  the  following  approximate  form  is  used: 

2 


y  =  62  (p')‘ 


(120) 


The  constant  6_  is  found  by  forcing  the  curve  (120)  to  intersect  the  dashed 

line  (118)  at  y  =  y  where  C  =0.9.  (The  C  was  chosen  by  trial  and  error 
o  o  o 

to  get  the  quadratic  y  in  (120)  as  close  as  possible  to  the  numerical  results 

within  the  operating  region.)  Namely  y  =  y  in  (118)  determines  the  abscissa 

1  o 

p'  of  the  intercept.  Finally  y  =  y^  and  the  p'  thus  determined  are  substituted 
in  (120)  to  yield 

6 2  =  1.125  (Bu(0)  )1/4.  (121) 

All  that  remains  is  to  write  the  wave  correlation  as  (55)  using  (116),  (120) ,  and 

(121) .  Equation  (55)  is  said  to  be  the  appropriate  one  if  p'  at  the  intercept 

is  greater  than  1/(k  L  ): 

m  o 

H  L  ^  3.3  (  B  (0)  )5/8.  (122) 

m  o  u 

If  (122)  is  not  satisfied,  then  the  form  (54)  is  to  be  used. 
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APPENDIX  III 

Relations  Between  n-dimensional  Power  Spectra 


This  section  provides  the  mathematical  relations  between  power 
spectral  densities  in  one,  two,  and  three  dimensions.  If  an  n-dimensional 
random  field  is  isotropic  and  homogeneous ,  then  a  one-dimensional  sample 
of  the  field  is  sufficient  to  give  the  power  spectrum  in  any  dimension  (s  n). 
Some  of  the  results  of  this  appendix  a^e  used  by  the  inverse  problem  of 
chapter  III  section  G. 

The  basis  for  the  relations  is  that  the  (homogeneous  but  not-  neces¬ 
sarily  isotropic)  covariance  B  (A)  is  invariant  with  dimension  in  that  it  is 
the  (1 , 1)  moment  of  the  joint  probability  density  W  (z^  a)  : 

B(a)  =  J  J  z1  z2  W(z1  ,z2;  A)  dzx  dz2  .  (123) 

-00 

However,  the  power  spectrum  (where  n  is  the  dimension)  is  the  n-th 


order  Fourier  transform  of  B  ( A)  : 

V(n)(K)  =,[../  B(  A)  e"ilC'^dZT  (124) 

—  OD 

where  ic  =  (k  ,K  ,  ...  ,  K  ) ,  A  =  (A,  .A,,  ...  ,  )*  and  dA  =  dA.  dA0  •  •  •  dA  . 

1  2  n  12  n  12  n 

If  isotropy  is  assumed: 

B(  A)  =  B(  o  )  (125) 

where  p  =  |  A  |  ,  and  the  Fourier  transform  (124)  may  be  reduced  to  a  one¬ 
dimensional  integral.  In  particular,  the  relations  for  n=  1,2,3  are: 

B  (  p )  =  —  f  cos(Kp)  V(1)(K)dK  (126) 

TT  0 
O 

=  f  J  (*  P)  V(2)  (K)  k  dK  (127) 

j  o 
o 

=  —  f  sin(Kp)  V(3)(K)KdK  (128) 

n  J 
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where  K  =  j  K  J  .  Since  the  three  transforms  possess  inverses,  they  provide 

a  means  of  going  between  power  spectral  densities  of  dimension  ns3. 

fnl 

In  some  cases,  more  direct  expressions  which  relate  V  [k)  may 
be  written  from  (126-  128).  First, 


V(3)(k)  =  -  -i- 


a  v 
a  k 


(i) 


(129) 


is  found  by  substituting  (128)  into  the  inverse  of  (126),  and  then  taking  the 
derivative  with  respect  to  K.  (Tatarski  [3]  found  the  same  result  by  an 
alternate  method.)  Second,  a  power  spectrum  in  a  particular  dimension  is 
related  to  the  spectrum  in  the  next  lower  dimension  by 


V(n-1)(K)  =  f  ■  ?  K  ■  V(n)(tc')d«'  (130) 

K  J  K'2  -  K4 

(n) 

(provided  is  Bessel- transformable).  Unfortunately  a  transform  inverse  to 

(130)  which  is  a  one-dimensional  integral  does  not  exist.  Hence  to  go  from 
to  V^n+1\  the  transforms  in  (126-128)  and  their  inverses  must  be  used. 
Those  transforms  are  in  fact  easier  to  use  than  (130)  since  Fourier  integral 
tables  are  available. 
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